Updated 5/7/25 by K. Lenderman - all code up to date & graphs generated - commented out outputs for winter sampling
This is an R Markdown Notebook. When you execute code within the notebook, the results appear beneath the code.
Need help? Try Stackoverflow: https://stackoverflow.com/tags/ggplot2
Attaching package: ‘dplyr’
The following objects are masked from ‘package:plyr’:
arrange, count, desc, failwith, id, mutate, rename, summarise, summarize
The following objects are masked from ‘package:stats’:
filter, lag
The following objects are masked from ‘package:base’:
intersect, setdiff, setequal, union
Loading required package: lattice
Attaching package: ‘lubridate’
The following objects are masked from ‘package:base’:
date, intersect, setdiff, union
Attaching package: ‘purrr’
The following object is masked from ‘package:plyr’:
compact
Attaching package: ‘tidyr’
The following object is masked from ‘package:reshape2’:
smiths
#data_all <- list.files(path = "Lab_Data_TissueProcessing/raw_data/Files_by_Month", # Identify all CSV files
# pattern = "*.csv", full.names = TRUE) %>%
#lapply(read_csv) %>% # Store all files in list
# bind_rows # Combine data sets into one data set
#data_all # Print data to RStudio console
#as.data.frame(data_all) # Convert tibble to data.frame
#Filtering NAs and unnecessary columns
#data_all <- data_all %>% filter(!is.na(date_collected))
#data_all <- select(data_all, -light_regime, -oyster_zone)
#write.csv(data_all, "Master_files/tissue_processing_all_data.csv", row.names=FALSE)
##### USE THIS CODE TO MERGE DATA FILES - some of the files have columns that are not the same format (dates specifically), which was causing issues in merging. This code below should solve that problem. It converts all date columns to dates and m/d/y format #####
file_paths <- list.files(path = "Lab_Data_TissueProcessing/raw_data/Files_by_Month",
pattern = "*.csv", full.names = TRUE)
data_all <- lapply(file_paths, function(file_path) {
read_csv(file_path) %>%
mutate(
date_collected = as.Date(date_collected,format = "%m/%d/%Y"),
date_processed = as.Date(date_processed,format = "%m/%d/%Y"),
date_davidsons = as.Date(date_davidsons,format = "%m/%d/%Y"),
date_etoh = as.Date(date_etoh,format = "%m/%d/%Y")
)
}) %>%
bind_rows
New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 38 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, oyster_zone, condition
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (20): light_regime, dissection_notes, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24...
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (18): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, ...27, ...28, .....
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 27── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (18): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17, ...18, ...19, ...20, ...21, ...22, .....
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 27── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (18): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17, ...18, ...19, ...20, ...21, ...22, .....
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (5): lab_id, site, oyster_zone, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (19): light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, .....
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (18): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, ...27, ...28, .....
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (20): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ......
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (20): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ......
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 27── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (18): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17, ...18, ...19, ...20, ...21, ...22, .....
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 35 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 46 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (5): lab_id, site, oyster_zone, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (19): light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, .....
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (19): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, date_etoh, ...21, ...22, ...23, ...24, ...25, ...26, ...27...
date (3): date_collected, date_processed, date_davidsons
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (20): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ......
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (20): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ......
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 27── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (17): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17, ...18, ...19, ...20, ...21, ...22, ...23, ...24, .....
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (5): lab_id, site, oyster_zone, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (19): light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, .....
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 998 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (19): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, date_etoh, ...21, ...22, ...23, ...24, ...25, ...26, ...27...
date (3): date_collected, date_processed, date_davidsons
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (21): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, date_etoh, ...21, ...22, ...23,...
date (3): date_collected, date_processed, date_davidsons
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (21): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, date_etoh, ...21, ...22, ...23,...
date (3): date_collected, date_processed, date_davidsons
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (5): lab_id, site, oyster_zone, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (19): light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, .....
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (19): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, date_etoh, ...21, ...22, ...23, ...24, ...25, ...26, ...27...
date (3): date_collected, date_processed, date_davidsons
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (21): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, date_etoh, ...21, ...22, ...23,...
date (3): date_collected, date_processed, date_davidsons
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (21): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, date_etoh, ...21, ...22, ...23,...
date (3): date_collected, date_processed, date_davidsons
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (5): lab_id, site, oyster_zone, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (19): light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, .....
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (18): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, ...27, ...28, .....
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (20): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ......
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (20): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ......
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (5): lab_id, site, oyster_zone, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (19): light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, .....
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (18): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, ...27, ...28, .....
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 36 Columns: 20── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (7): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (20): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ......
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (5): lab_id, site, oyster_zone, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (19): light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, .....
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (18): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, ...27, ...28, .....
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 27── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (17): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17, ...18, ...19, ...20, ...21, ...22, ...23, ...24, .....
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
as.data.frame(data_all) # Convert tibble to data.frame
data_all <- data_all %>% filter(!is.na(date_collected))
data_all <- select(data_all, -light_regime, -oyster_zone)
View(data_all)
write.csv(data_all, "Master_files/tissue_processing_all_data.csv", row.names=FALSE)
#A wrong value caliper input was identified for the height of sample 0923GOLD_23. This code is removing that value from the data set as we cannot conclude what this original value was. The value is 8.62. This will cause this individual to fall out of the dataset when standardized to length. This code does not completely remove the individual from the dataset.
data_all$height_mm[data_all$height_mm == "0"] <- NA
data_all <- data_all %>% dplyr::mutate(date_collected= as.Date(date_collected), month = month(date_collected))
data_all <- data_all %>% dplyr::mutate(date_collected= as.Date(date_collected), year = year(date_collected))
#changing numeric month to month name
data_all$month <- factor(data_all$month, levels = c("1","2","3","4","5","6","7", "8", "9", "10", "11", "12"),
labels=c("Jan","Feb", "March","April","May", "June", "July", "Aug", "Sept", "Oct", "Nov", "Dec"))
data_all
data_all<- data_all %>%
dplyr::mutate(condition_score = recode(condition, "1_very_good" = 1, "2_good" = 2, "3_good_minus"= 3, "4_fair_plus"= 4, "5_fair"= 5,"6_fair_minus"=6,"7_poor_plus"=7, "8_poor"= 8, "9_very_poor"= 9))
Warning: There was 1 warning in `dplyr::mutate()`.
ℹ In argument: `condition_score = recode(...)`.
Caused by warning:
! Unreplaced values treated as NA as `.x` is not compatible.
Please specify replacements exhaustively or supply `.default`.
head(data_all)
NA
NA
#This chunk of code is removing 0723LAUR_20 and 0723LAUR_26 from the datasheet as they have been identified as spat on shell to avoid bias in the data. During this sample collection there were animals that were significantly smaller than the single set oysters. These individuals should be removed from all monthly sampling related datasheets including disease analysis. All tissue amples will be disgarded.
data_all <- data_all %>%
subset(lab_id != "0723LAUR_20") %>%
subset(lab_id != "0723LAUR_26")
data_all
NA
st_height <- summarySE(data_all%>% filter(!is.na(height_mm)), measurevar="height_mm", groupvars=c("site", "date_collected"))
st_height
#Calculate completeness for QC
st_height$Completeness <- st_height$N /30
st_height
write.csv(st_height, "Lab_Data_TissueProcessing/output\\Completeness_tissue_processing_data.csv", row.names=FALSE)
#Mean Height
ggplot(data=data_all, aes(x=site, y=height_mm, fill=site)) +
geom_boxplot()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Mean Shell Height ", x ="site", y = "Mean Shell Height (mm)") + facet_wrap(.~year)
#Body condition
ggplot(data=data_all, aes(x=site, y=condition_score, fill=site)) +
geom_boxplot()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Mean Condition Score ", x ="site", y = "Mean Body Condition Score")+ scale_y_reverse()+ facet_wrap(.~year)
mean_body_condition <- data_all %>%
dplyr::group_by(site, month, year)%>%
dplyr::summarize(mean_bsc = mean(condition_score, na.rm = TRUE))
`summarise()` has grouped output by 'site', 'month'. You can override using the `.groups` argument.
mean_body_condition
#for overlay graph
#mean_bcs_ashc <- mean_body_condition %>%
#filter(site == "ASHC") %>%
#filter(year == "2024")
#Mean body condtion scores - 2023
data_2023 <- data_all%>% filter(year=="2023")
ggplot(data=data_2023, aes(x=month, y=condition_score, group = month, fill = site)) +
geom_boxplot()+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = -35, vjust = 0.5, hjust=1))+
labs(title="Mean Condition Score - 2023", x ="month", y = "Mean Body Condition Score")+ scale_y_reverse()+facet_wrap(~site)
#Mean body condtion scores - 2024
data_2024 <- data_all%>% filter(year=="2024")
ggplot(data=data_2024, aes(x=month, y=condition_score, group = month, fill = site)) +
geom_boxplot()+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = -35, vjust = 0.5, hjust=1))+
labs(title="Mean Condition Score - 2024", x ="month", y = "Mean Body Condition Score")+ scale_y_reverse()+facet_wrap(~site)
#Mean body condition scores - 2025
data_2025 <- data_all%>% filter(year=="2025")
ggplot(data=data_2025, aes(x=month, y=condition_score, group = month, fill = site)) +
geom_boxplot()+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = -35, vjust = 0.5, hjust=1))+
labs(title="Mean Condition Score - 2025", x ="month", y = "Mean Body Condition Score")+ scale_y_reverse()+facet_wrap(~site)
#Proportions graph Body condition scores
#Proportions graph Body condition scores - 2023
df_BCS_proportions2023<- data_2023 %>%
mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ "1", condition_score == 2 ~"2",condition_score == 3 ~ "3", condition_score == 4 ~"4",condition_score == 5 ~ "5", condition_score == 6 ~"6",condition_score == 7 ~ "7", condition_score == 8 ~"8", condition_score ==9 ~"9", TRUE ~ as.character(condition_score))) %>%
group_by(site, month, year, condition_score_bin) %>% dplyr ::summarise(Count= n()) %>%
ungroup() %>%
mutate(Proportion = Count/sum(Count))
`summarise()` has grouped output by 'site', 'month', 'year'. You can override using the `.groups` argument.
df_BCS_proportions2023<- na.omit(df_BCS_proportions2023)
BCS_proportion_all_2023<- ggplot(data=df_BCS_proportions2023, aes(x=month, y= Proportion, fill=condition_score_bin)) +
geom_bar(width = .5, stat="identity", position = "fill", colour = "black")+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
labs(title="Proportion of body condition scores - 2023", x ="month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
theme(axis.text=element_text(size=12))+
scale_fill_brewer(palette = "Blues", direction = -1)+
facet_wrap(~ site)
BCS_proportion_all_2023
#pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/all_BCS_proportion_2023.pdf"), height = 7, width = 13)
#print(BCS_proportion_all_2023)
#dev.off()
#Proportions graph Body condition scores - 2024
df_BCS_proportions2024<- data_2024 %>%
mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ "1", condition_score == 2 ~"2",condition_score == 3 ~ "3", condition_score == 4 ~"4",condition_score == 5 ~ "5", condition_score == 6 ~"6",condition_score == 7 ~ "7", condition_score == 8 ~"8", condition_score ==9 ~"9", TRUE ~ as.character(condition_score))) %>%
group_by(site, month, year, condition_score_bin) %>% dplyr ::summarise(Count= n()) %>%
ungroup() %>%
mutate(Proportion = Count/sum(Count))
`summarise()` has grouped output by 'site', 'month', 'year'. You can override using the `.groups` argument.
df_BCS_proportions2024<- na.omit(df_BCS_proportions2024)
BCS_proportion_all_2024<- ggplot(data=df_BCS_proportions2024, aes(x=month, y= Proportion, fill=condition_score_bin)) +
geom_bar(width = .5, stat="identity", position = "fill", colour = "black")+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
labs(title="Proportion of body condition scores - 2024", x ="month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
theme(axis.text=element_text(size=12))+
scale_fill_brewer(palette = "Blues", direction = -1)+
facet_wrap(~ site)
BCS_proportion_all_2024
#pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/all_BCS_proportion_2024.pdf"), height = 7, width = 13)
#print(BCS_proportion_all_2024)
#dev.off()
#Proportions graph Body condition scores - 2025
df_BCS_proportions2025<- data_2025 %>%
mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ "1", condition_score == 2 ~"2",condition_score == 3 ~ "3", condition_score == 4 ~"4",condition_score == 5 ~ "5", condition_score == 6 ~"6",condition_score == 7 ~ "7", condition_score == 8 ~"8", condition_score ==9 ~"9", TRUE ~ as.character(condition_score))) %>%
group_by(site, month, year, condition_score_bin) %>% dplyr ::summarise(Count= n()) %>%
ungroup() %>%
mutate(Proportion = Count/sum(Count))
`summarise()` has grouped output by 'site', 'month', 'year'. You can override using the `.groups` argument.
df_BCS_proportions2025<- na.omit(df_BCS_proportions2025)
BCS_proportion_all_2025<- ggplot(data=df_BCS_proportions2025, aes(x=month, y= Proportion, fill=condition_score_bin)) +
geom_bar(width = .5, stat="identity", position = "fill", colour = "black")+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
labs(title="Proportion of body condition scores - 2025", x ="month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
theme(axis.text=element_text(size=12))+
scale_fill_brewer(palette = "Blues", direction = -1)+
facet_wrap(~ site)#, scales = "free")
BCS_proportion_all_2025
pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/all_BCS_proportion_2025.pdf"), height = 7, width = 13)
print(BCS_proportion_all_2025)
dev.off()
png
2
#Winter sampling
# winter_bcs <- rbind(df_BCS_proportions2024, df_BCS_proportions2025)
#
# winter_bcs <- winter_bcs %>% filter(month == "Dec"|month == "Jan"|month == "Feb")
#
# BCS_proportion_winter<- ggplot(data=winter_bcs, aes(x=factor (month, level=c('Dec', 'Jan', 'Feb')), y= Proportion, fill=condition_score_bin)) +
# geom_bar(width = .5, stat="identity", position = "fill", colour = "black")+
# theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
# theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
# labs(title="Proportion of body condition scores - Winter Sampling 2024-2025", x ="Month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
# theme(axis.text=element_text(size=12))+
# scale_fill_brewer(palette = "Blues", direction = -1)+
# facet_wrap(~ site)
#
# BCS_proportion_winter
# Proportion per month of intensity scores at three sites during 2023 & 2024 & 2025
df_bcs_proportions_new <- rbind(df_BCS_proportions2023, df_BCS_proportions2024, df_BCS_proportions2025)
df_bcs_proportions_new<- df_bcs_proportions_new %>% filter(!site=="LAUR")
ggplot(data=df_bcs_proportions_new, aes(x=month, y= Proportion, fill=condition_score_bin)) +
geom_bar(width = .5, stat="identity", position = "fill",colour = "black")+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
labs(fill ="Body condition scores")+
labs(title="Proportion of body condition scores", x ="month", y = "Proportion")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
theme(axis.text=element_text(size=12))+
scale_fill_brewer(palette = "Blues", direction = -1)+
#scale_x_continuous("Month", breaks = c(3,4,5,6,7,8,9,10))+
facet_grid(year ~site)
#Shell Pathology - all data #Includes 2023 data and all data that is in the 2024 format
shell_path <- list.files(path = "Lab_Data_TissueProcessing/raw_data/shell_pathology", # Identify all CSV files
pattern = "*.csv", full.names = TRUE) %>%
lapply(read_csv) %>% # Store all files in list
bind_rows # Combine data sets into one data set
Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 26── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (2): lab_id, pathology_notes
lgl (24): boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discoloration, horn_add, ...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 26── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (2): lab_id, pathology_notes
lgl (24): boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discoloration, horn_add, ...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 26── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (25): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 26── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (25): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 26── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (25): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 26── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (25): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 45 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (2): lab_id, pathology_notes
lgl (12): boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discoloration, horn_add, ...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 15── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (14): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 28 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 30 Columns: 15── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (2): lab_id, ...15
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 26── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): lab_id
lgl (25): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 838 Columns: 11── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (2): lab_id, pathology_notes
lgl (9): boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, horn_add, shell_scarring, cyst_abscess
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
shell_path <- shell_path %>% filter(!is.na(lab_id))
shell_path <- shell_path %>% select(-c("...14":"...26"))
write.csv(shell_path, "Master_files/shell_pathology_all_data.csv", row.names=FALSE)
#Writing summary for shell pathology
shell_path_summary <- shell_path%>% separate(lab_id, c('Date_Site', 'ID'))
shell_path_summary <- select(shell_path_summary, -c("pathology_notes","ID",))
shell_path_summary$boring_sponge <- as.integer(as.logical(shell_path_summary$boring_sponge))
shell_path_summary$polydora <- as.integer(as.logical(shell_path_summary$polydora))
shell_path_summary$conchiolin_mod_severe <- as.integer(as.logical(shell_path_summary$conchiolin_mod_severe))
shell_path_summary$mud_blister <- as.integer(as.logical(shell_path_summary$mud_blister))
shell_path_summary$pea_crab <- as.integer(as.logical(shell_path_summary$pea_crab))
shell_path_summary$gill_erosion <- as.integer(as.logical(shell_path_summary$gill_erosion))
shell_path_summary$pale_digestive <- as.integer(as.logical(shell_path_summary$pale_digestive))
shell_path_summary$discoloration <- as.integer(as.logical(shell_path_summary$discoloration))
shell_path_summary$horn_add <- as.integer(as.logical(shell_path_summary$horn_add))
shell_path_summary$cyst_abscess <- as.integer(as.logical(shell_path_summary$cyst_abscess))
shell_path_summary$tumor <- as.integer(as.logical(shell_path_summary$tumor))
shell_path_summary$oyster_drill <- as.integer(as.logical(shell_path_summary$oyster_drill))
shell_path_summary$boring_snail <- as.integer(as.logical(shell_path_summary$boring_snail))
shell_path_summary$shell_scarring <- as.integer(as.logical(shell_path_summary$shell_scarring))
shell_path_count<- shell_path_summary %>%
dplyr::group_by(Date_Site) %>%
dplyr::summarize(Sample_count = n(),
Boring_sponge =sum(boring_sponge),
Polydora = sum(polydora),
Conchiolin =sum(conchiolin_mod_severe),
Mud_blister =sum(mud_blister),
Pea_crab =sum(pea_crab),
Gill_erosion =sum(gill_erosion),
Pale_digestive =sum(pale_digestive),
Discoloration =sum(discoloration),
Horn =sum(horn_add),
Cyst =sum(cyst_abscess),
Tumor =sum(tumor),
Oyster_drill =sum(oyster_drill),
Boring_snail =sum(boring_snail),
Shell_scarring =sum(shell_scarring),
Pathogen_count = sum(boring_sponge,polydora,conchiolin_mod_severe, mud_blister, pea_crab,
gill_erosion, pale_digestive, discoloration, horn_add, cyst_abscess,
tumor, oyster_drill, boring_snail, shell_scarring, na.rm = TRUE)) %>%
ungroup()
write.csv(shell_path_count, "Lab_Data_TissueProcessing/output/shell_pathology_counts.csv", row.names=FALSE)
shell_path_summary
df_ASHC<- data_all%>%
filter(site=="ASHC")
df_ASHC
## Shell Height ##
st_height_ASHC <- summarySE(df_ASHC, measurevar="height_mm", groupvars=c("date_collected"))
st_height_ASHC
## Body Condition ##
#Excludes April and May due to scoring change. These months are scored categorically 'fat, medium, watery'.
st_condition_ASHC<- summarySE(df_ASHC, measurevar = "condition_score", groupvars = c("date_collected"))
st_condition_ASHC
mean_condition_ASHC <- st_condition_ASHC %>%
filter(!is.na(condition_score)) %>% # Filter out rows where condition_score is NA
summarize(mean_bcs = mean(condition_score, na.rm = TRUE))
mean_condition_ASHC
NA
ggplot(data=df_ASHC, aes(x=month, y=height_mm, group= month)) +
geom_boxplot()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title=" Ash Creek Mean Shell Height ", x ="month", y = "Mean Shell Height (mm)")+facet_wrap(.~year)
#Condition distribution across all sample months - 2023
df_ASHC_2023 <- df_ASHC%>% filter(!year=="2024")
ASHC_BCS_dist2023 <-ggplot(data=df_ASHC_2023, aes(x= condition)) +
geom_bar()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Ash Creek Body Condition Score - 2023", x ="Condition categorization")+
facet_wrap(~ month, scales = "free")
ASHC_BCS_dist2023
#pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/ASHC_BCS_dist_2023.pdf"), height = 7, width = 13)
#print(ASHC_BCS_dist2023)
#dev.off()
#Condition distribution across all sample months - 2024
df_ASHC_2024 <- df_ASHC%>% filter(!year=="2023")
ASHC_BCS_dist2024 <-ggplot(data=df_ASHC_2024, aes(x= condition)) +
geom_bar()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Ash Creek Body Condition Score - 2024", x ="Condition categorization")+
facet_wrap(~ month, scales = "free")
ASHC_BCS_dist2024
#pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/ASHC_BCS_dist_2024.pdf"), height = 7, width = 13)
#print(ASHC_BCS_dist2024)
#dev.off()
#Condition distribution across all sample months - 2025
df_ASHC_2025 <- df_ASHC%>% filter(year=="2025")
ASHC_BCS_dist2025 <-ggplot(data=df_ASHC_2025, aes(x= condition)) +
geom_bar()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Ash Creek Body Condition Score - 2025", x ="Condition categorization")+
facet_wrap(~ month, scales = "free")
ASHC_BCS_dist2025
pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/ASHC_BCS_dist_2025.pdf"), height = 7, width = 13)
print(ASHC_BCS_dist2025)
dev.off()
png
2
#Mean Body condition per month
ggplot(data=df_ASHC, aes(x= month, y= condition_score, group = month)) +
geom_boxplot()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Ash Creek Mean Body Condition Score ", x ="month", y= " condition score (1-9)") + scale_y_reverse()+facet_wrap(.~year)
#ASHC Proportions graph Body condition scores - all years
ASHC_BCS_proportions<- df_ASHC %>%
mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ "1", condition_score == 2 ~"2",condition_score == 3 ~ "3", condition_score == 4 ~"4",condition_score == 5 ~ "5", condition_score == 6 ~"6",condition_score == 7 ~ "7", condition_score == 8 ~"8", condition_score ==9 ~"9", TRUE ~ as.character(condition_score))) %>%
group_by(site, month, condition_score_bin, year) %>% dplyr::summarise(Count= n()) %>%
ungroup() %>%
mutate(Proportion = Count/sum(Count))
`summarise()` has grouped output by 'site', 'month', 'condition_score_bin'. You can override using the `.groups` argument.
ASHC_BCS_proportions<- na.omit(ASHC_BCS_proportions)
ASHC_BCS_proportions
BCS_proportion_ASHC<- ggplot(data=ASHC_BCS_proportions, aes(x=month, y= Proportion, fill=condition_score_bin)) +
geom_bar(width = .5, stat="identity", position = "fill")+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
labs(title="Ash Creek Proportion of body condition scores June- November", x ="month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
theme(axis.text=element_text(size=12))+
facet_wrap(~ year)
#scale_fill_brewer() +
#facet_wrap(~ site)
BCS_proportion_ASHC
pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/ASHC_BCS_proportion.pdf"), height = 7, width = 13)
print(BCS_proportion_ASHC)
dev.off()
png
2
df_ASHC
ASHC_precent_greater_3 <- df_ASHC %>%
dplyr::group_by(month, site, year) %>%
dplyr::summarise(Percentage = mean(condition_score <= 3)*100)
`summarise()` has grouped output by 'month', 'site'. You can override using the `.groups` argument.
ASHC_precent_greater_3
ASHC_BCS_percentage <- ASHC_precent_greater_3%>%
#filter(year =="2024")%>%
ggplot(aes(x = month, y = Percentage)) +
#geom_bar(width = 0.5, stat = "identity", position = "fill") +
geom_col()+
theme_bw() +
theme(
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1),
axis.title.y = element_text(size = rel(1.3), angle = 90),
axis.title.x = element_text(size = rel(1.3), angle = 0),
axis.text = element_text(size = 12)
) +
labs(
title = "Ash Creek % of body condition scores >= 3",
x = "month",
y = "Percentage of body condition scores >= 3"
) +
facet_wrap(~year)
#scale_x_continuous( breaks = seq(5,12, by =1) )+
#ylim(0,60)+
#scale_y_continuous(limits = c(0,100), breaks = seq(0,100, by = 10))
ASHC_BCS_percentage
pdf(paste0(path = "Lab_Data_TissueProcessing/output", "/ASHC_BCS_percentage.pdf"),height = 7, width = 13)
print(ASHC_BCS_percentage)
dev.off()
png
2
df_FENC<- data_all%>%
filter(site=="FENC")
df_FENC
## Shell Height ##
st_height_FENC <- summarySE(df_FENC, measurevar="height_mm", groupvars=c("date_collected"))
st_height_FENC
## Body Condition ##
#Excludes April and May due to scoring change. These months are scored categorically 'fat, medium, watery'.
st_condition_FENC<- summarySE(df_FENC%>% filter(!is.na(condition_score)), measurevar = "condition_score", groupvars = c("date_collected"))
st_condition_FENC
NA
ggplot(data=df_FENC, aes(x=month, y=height_mm, group=month)) +
geom_boxplot()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title=" Fence Creek Mean Shell Height ", x ="month", y = "Mean Shell Height (mm)")+ facet_wrap(.~year)
#Condition distribution across all sample months - 2023
df_FENC.2023 <- df_FENC%>% filter(!year=="2024")
ggplot(data=df_FENC.2023, aes(x= condition)) +
geom_bar()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Fence Creek Body Condition index - 2023", x ="Condition categorization")+
facet_wrap(~ month, scales = "free")
#Condition distribution across all sample months - 2024
df_FENC.2024 <- df_FENC%>% filter(!year=="2023")
ggplot(data=df_FENC.2024, aes(x= condition)) +
geom_bar()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Fence Creek Body Condition index - 2024", x ="Condition categorization")+
facet_wrap(~ month, scales = "free")
#Mean Body condition per month
ggplot(data=df_FENC, aes(x= month, y= condition_score, group = month)) +
geom_boxplot()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Fence Creek Mean Body Condition Score", x ="month", y= " condition score (1-9)")+ scale_y_reverse() + facet_wrap(~ year, scales = "free")
#FENC Proportions graph Body condition scores - all years
FENC_BCS_proportions<- df_FENC %>%
mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ "1", condition_score == 2 ~"2",condition_score == 3 ~ "3", condition_score == 4 ~"4",condition_score == 5 ~ "5", condition_score == 6 ~"6",condition_score == 7 ~ "7", condition_score == 8 ~"8", condition_score ==9 ~"9", TRUE ~ as.character(condition_score))) %>%
group_by(site, month, condition_score_bin, year) %>% dplyr::summarise(Count= n()) %>%
ungroup() %>%
mutate(Proportion = Count/sum(Count))
`summarise()` has grouped output by 'site', 'month', 'condition_score_bin'. You can override using the `.groups` argument.
FENC_BCS_proportions<- na.omit(FENC_BCS_proportions)
FENC_BCS_proportions
BCS_proportion_FENC<- ggplot(data=FENC_BCS_proportions, aes(x=month, y= Proportion, fill=condition_score_bin)) +
geom_bar(width = .5, stat="identity", position = "fill")+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
labs(title="Fence Creek Proportion of body condition scores", x ="month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
theme(axis.text=element_text(size=12))+
facet_wrap(~ year)
#scale_fill_brewer() +
#facet_wrap(~ site)
BCS_proportion_FENC
pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/FENC_BCS_proportion.pdf"), height = 7, width = 13)
print(BCS_proportion_FENC)
dev.off()
png
2
#removing 0823FENC_28 bcs was NA
df_FENC <- df_FENC %>% drop_na(condition)
FENC_precent_greater_3 <- df_FENC %>%
dplyr::group_by(month, year) %>%
dplyr::summarise(Percentage = mean(condition_score <= 3)*100)
`summarise()` has grouped output by 'month'. You can override using the `.groups` argument.
FENC_precent_greater_3
FENC_BCS_percentage <- ggplot(FENC_precent_greater_3,aes(x = month, y = Percentage)) +
#geom_bar(width = 0.5, stat = "identity", position = "fill") +
geom_col()+
theme_bw() +
theme(
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1),
axis.title.y = element_text(size = rel(1.3), angle = 90),
axis.title.x = element_text(size = rel(1.3), angle = 0),
axis.text = element_text(size = 12)
) +
labs(
title = "Fence Creek % of body condition scores >= 3",
x = "month",
y = "Percentage of body condition scores >= 3"
) +
facet_wrap(~year)
#scale_x_continuous( breaks = seq(5,12, by =1) )+
#ylim(0,60)+
#scale_y_continuous(limits = c(0,100), breaks = seq(0,100, by = 10))
FENC_BCS_percentage
pdf(paste0(path = "Lab_Data_TissueProcessing/output", "/FENC_BCS_percentage.pdf"),height = 7, width = 13)
print(FENC_BCS_percentage)
dev.off()
png
2
df_GOLD<- data_all%>%
filter(site=="GOLD")
df_GOLD <- df_GOLD %>% filter(!is.na(height_mm))
#2023 data
df_GOLD2023 <- df_GOLD%>% filter(!year=="2024")
#2024 data
df_GOLD2024 <- df_GOLD%>% filter(!year=="2023")
df_GOLD$height_mm<-as.numeric(df_GOLD$height_mm)
mean_shell_height <- df_GOLD %>% mutate(year = year(date_collected)) %>% group_by(year)%>% dplyr::summarise(mean_height = mean(height_mm)) #summarySE(measurevar="height_mm", groupvars=c("year"))
mean_shell_height
## Shell Height ##
st_height_GOLD <- summarySE(df_GOLD, measurevar="height_mm", groupvars=c("date_collected"))
st_height_GOLD
## Body Condition ##
st_condition_GOLD<- summarySE(df_GOLD, measurevar = "condition_score", groupvars = c("date_collected"))
st_condition_GOLD
NA
#GOLD Proportions graph Body condition scores
GOLD_BCS_proportions<- df_GOLD %>%
mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ "1", condition_score == 2 ~"2",condition_score == 3 ~ "3", condition_score == 4 ~"4",condition_score == 5 ~ "5", condition_score == 6 ~"6",condition_score == 7 ~ "7", condition_score == 8 ~"8", condition_score ==9 ~"9", TRUE ~ as.character(condition_score))) %>%
group_by(site, month, year, condition_score_bin) %>% dplyr::summarise(Count= n()) %>%
ungroup() %>%
mutate(Proportion = Count/sum(Count))
`summarise()` has grouped output by 'site', 'month', 'year'. You can override using the `.groups` argument.
GOLD_BCS_proportions<- na.omit(GOLD_BCS_proportions)
BCS_proportion_GOLD <- ggplot(data=GOLD_BCS_proportions, aes(x=month, y= Proportion, fill=condition_score_bin)) +
geom_bar(width = .5, stat="identity", position = "fill", colour = "black")+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
labs(title="Gold Star Beach Proportion of body condition scores", x ="month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
theme(axis.text=element_text(size=12))+
scale_fill_brewer(palette = "Blues", direction = -1)+
facet_wrap(~ year)
BCS_proportion_GOLD
pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/GOLD_BCS_proportion.pdf"), height = 7, width = 13)
print(BCS_proportion_GOLD)
dev.off()
png
2
df_GOLD
GOLD_precent_greater_3 <- df_GOLD %>%
dplyr::group_by(month, site, year) %>%
dplyr::summarise(Percentage = mean(condition_score <= 3)*100)
`summarise()` has grouped output by 'month', 'site'. You can override using the `.groups` argument.
GOLD_precent_greater_3
GOLD_BCS_percentage <- GOLD_precent_greater_3 %>%
ggplot(aes(x = month, y = Percentage)) +
#geom_bar(width = 0.5, stat = "identity", position = "fill") +
geom_col()+
theme_bw() +
theme(
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1),
axis.title.y = element_text(size = rel(1.3), angle = 90),
axis.title.x = element_text(size = rel(1.3), angle = 0),
axis.text = element_text(size = 12)
) +
labs(
title = "Gold Star Beach % of body condition scores >= 3",
x = "month",
y = "Percentage of body condition scores >= 3"
) +facet_wrap(~year)
#+ scale_x_continuous( breaks = seq(5,12, by =1) )+
#ylim(0,60)+
#scale_y_continuous(limits = c(0,100), breaks = seq(0,100, by = 10))
GOLD_BCS_percentage
pdf(paste0(path = "Lab_Data_TissueProcessing/output", "/GOLD_BCS_percentage.pdf"),height = 7, width = 13)
print(GOLD_BCS_percentage)
dev.off()
png
2
ggplot(data=df_GOLD, aes(x=month, y=height_mm, group= month)) +
geom_boxplot()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Goldstar Beach Mean Shell Height ", x ="month", y = "Mean Shell Height (mm)") + #scale_x_continuous(limits= c(4,11), breaks = seq(5,10, by =1))+
geom_smooth(method = "lm", se = FALSE)+ facet_wrap(~year)
#Condition distribution across all sample months - 2023
ggplot(data=df_GOLD2023, aes(x= condition)) +
geom_bar()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Goldstar Beach Body Condition Score - 2023 ", x ="Condition categorization")+
facet_wrap(~ month, scales = "free")
#Condition distribution across all sample months - 2024
ggplot(data=df_GOLD2024, aes(x= condition)) +
geom_bar()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Goldstar Beach Body Condition Score - 2024", x ="Condition categorization")+
facet_wrap(~ month, scales = "free")
#Mean Body condition per month
ggplot(data=df_GOLD, aes(x= month, y= condition_score, group = month)) +
geom_boxplot()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
labs(title="Goldstar Beach Mean Body Condition Score ", x ="month", y= " condition score (1-9)")+ scale_y_reverse() + facet_wrap(~year)
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
<!-- rnb-source-begin eyJkYXRhIjoiYGBgclxuZGZfTEFVUjwtIGRhdGFfYWxsJT4lXG4gIGZpbHRlcihzaXRlPT1cIkxBVVJcIilcblxuIyMgU2hlbGwgSGVpZ2h0ICMjXG5zdF9oZWlnaHRfTEFVUiA8LSBzdW1tYXJ5U0UoZGZfTEFVUiwgbWVhc3VyZXZhcj1cImhlaWdodF9tbVwiLCBncm91cHZhcnM9YyhcImRhdGVfY29sbGVjdGVkXCIpKVxuc3RfaGVpZ2h0X0xBVVIgXG5cbiMjIEJvZHkgQ29uZGl0aW9uICMjXG4jIyBFeGNsdWRlcyBNYXkgZHVlIHRvIHNjb3Jpbmcgc3lzdGVtIGNoYW5nZS4gTWF5IHNjb3JlZCBmYXQsIG1lZGl1bSwgd2F0ZXJ5LiBcbnN0X2NvbmRpdGlvbl9MQVVSPC0gc3VtbWFyeVNFKGRmX0xBVVIsIG1lYXN1cmV2YXIgPSBcImNvbmRpdGlvbl9zY29yZVwiLCBncm91cHZhcnMgPSBjKFwiZGF0ZV9jb2xsZWN0ZWRcIikpXG5zdF9jb25kaXRpb25fTEFVUlxuYGBgIn0= -->
```r
df_LAUR<- data_all%>%
filter(site==\LAUR\)
## Shell Height ##
st_height_LAUR <- summarySE(df_LAUR, measurevar=\height_mm\, groupvars=c(\date_collected\))
st_height_LAUR
## Body Condition ##
## Excludes May due to scoring system change. May scored fat, medium, watery.
st_condition_LAUR<- summarySE(df_LAUR, measurevar = \condition_score\, groupvars = c(\date_collected\))
st_condition_LAUR
```
<!-- rnb-source-end -->
```r
df_LAUR<- data_all%>%
filter(site==\LAUR\)
## Shell Height ##
st_height_LAUR <- summarySE(df_LAUR, measurevar=\height_mm\, groupvars=c(\date_collected\))
st_height_LAUR
## Body Condition ##
## Excludes May due to scoring system change. May scored fat, medium, watery.
st_condition_LAUR<- summarySE(df_LAUR, measurevar = \condition_score\, groupvars = c(\date_collected\))
st_condition_LAUR
<!-- rnb-source-end -->
<!-- rnb-output-end -->
<!-- rnb-chunk-end -->
<!-- rnb-text-begin -->
<!-- rnb-text-end -->
<!-- rnb-chunk-begin -->
<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVaMmR3Ykc5MEtHUmhkR0U5WkdaZlRFRlZVaXdnWVdWektIZzliVzl1ZEdnc0lIazlhR1ZwWjJoMFgyMXRMQ0JuY205MWNEMGdiVzl1ZEdncEtTQXJYRzRnSUdkbGIyMWZZbTk0Y0d4dmRDZ3BLeUFnSTNOallXeGxYMlpwYkd4ZmJXRnVkV0ZzS0haaGJIVmxjejFqS0Z3aVptOXlaWE4wWjNKbFpXNWNJaXhjSW05eVlXNW5aVndpTENCY0luQjFjbkJzWlZ3aUtTa3JYRzRnSUhSb1pXMWxYMkozS0NrZ0t5QWdkR2hsYldVb2NHRnVaV3d1WjNKcFpDNXRZV3B2Y2lBOUlHVnNaVzFsYm5SZllteGhibXNvS1N3Z2NHRnVaV3d1WjNKcFpDNXRhVzV2Y2lBOUlHVnNaVzFsYm5SZllteGhibXNvS1NrcklGeHVJQ0FqZEdobGJXVW9ZWGhwY3k1MFpYaDBMbmdnUFNCbGJHVnRaVzUwWDNSbGVIUW9ZVzVuYkdVZ1BTQTVNQ3dnZG1wMWMzUWdQU0F3TGpVc0lHaHFkWE4wUFRFcEtTdGNiaUFnYkdGaWN5aDBhWFJzWlQxY0lpQk1ZWFZ5Wld3Z1NHOXNiRzkzSUUxbFlXNGdVMmhsYkd3Z1NHVnBaMmgwSUMwZ01qQXlNMXdpTENCNElEMWNJbTF2Ym5Sb1hDSXNJSGtnUFNCY0lrMWxZVzRnVTJobGJHd2dTR1ZwWjJoMElDaHRiU2xjSWlsY2JseHVYRzVnWUdBaWZRPT0gLS0+XG5cbmBgYHJcbmdncGxvdChkYXRhPWRmX0xBVVIsIGFlcyh4PW1vbnRoLCB5PWhlaWdodF9tbSwgZ3JvdXA9IG1vbnRoKSkgK1xuICBnZW9tX2JveHBsb3QoKSsgICNzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXM9YyhcXGZvcmVzdGdyZWVuXFwsXFxvcmFuZ2VcXCwgXFxwdXJwbGVcXCkpK1xuICB0aGVtZV9idygpICsgIHRoZW1lKHBhbmVsLmdyaWQubWFqb3IgPSBlbGVtZW50X2JsYW5rKCksIHBhbmVsLmdyaWQubWlub3IgPSBlbGVtZW50X2JsYW5rKCkpKyBcbiAgI3RoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gOTAsIHZqdXN0ID0gMC41LCBoanVzdD0xKSkrXG4gIGxhYnModGl0bGU9XFwgTGF1cmVsIEhvbGxvdyBNZWFuIFNoZWxsIEhlaWdodCAtIDIwMjNcXCwgeCA9XFxtb250aFxcLCB5ID0gXFxNZWFuIFNoZWxsIEhlaWdodCAobW0pXFwpXG5cblxuYGBgXG5cbjwhLS0gcm5iLXNvdXJjZS1lbmQgLS0+XG4ifQ== -->
ggplot(data=df_LAUR, aes(x=month, y=height_mm, group= month)) +
geom_boxplot()+ #scale_fill_manual(values=c(\forestgreen\,\orange\, \purple\))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title=\ Laurel Hollow Mean Shell Height - 2023\, x =\month\, y = \Mean Shell Height (mm)\)
<!-- rnb-output-end -->
<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVZR0JnY2x4dVoyZHdiRzkwS0dSaGRHRTlaR1pmVEVGVlVpd2dZV1Z6S0hnOWJXOXVkR2dzSUhrOWFHVnBaMmgwWDIxdExDQm5jbTkxY0QwZ2JXOXVkR2dwS1NBclhHNGdJR2RsYjIxZlltOTRjR3h2ZENncEt5QWdJM05qWVd4bFgyWnBiR3hmYldGdWRXRnNLSFpoYkhWbGN6MWpLRnhjWm05eVpYTjBaM0psWlc1Y1hDeGNYRzl5WVc1blpWeGNMQ0JjWEhCMWNuQnNaVnhjS1NrclhHNGdJSFJvWlcxbFgySjNLQ2tnS3lBZ2RHaGxiV1VvY0dGdVpXd3VaM0pwWkM1dFlXcHZjaUE5SUdWc1pXMWxiblJmWW14aGJtc29LU3dnY0dGdVpXd3VaM0pwWkM1dGFXNXZjaUE5SUdWc1pXMWxiblJmWW14aGJtc29LU2tySUZ4dUlDQWpkR2hsYldVb1lYaHBjeTUwWlhoMExuZ2dQU0JsYkdWdFpXNTBYM1JsZUhRb1lXNW5iR1VnUFNBNU1Dd2dkbXAxYzNRZ1BTQXdMalVzSUdocWRYTjBQVEVwS1N0Y2JpQWdiR0ZpY3loMGFYUnNaVDFjWENCTVlYVnlaV3dnU0c5c2JHOTNJRTFsWVc0Z1UyaGxiR3dnU0dWcFoyaDBJQzBnTWpBeU0xeGNMQ0I0SUQxY1hHMXZiblJvWEZ3c0lIa2dQU0JjWEUxbFlXNGdVMmhsYkd3Z1NHVnBaMmgwSUNodGJTbGNYQ2xjYmx4dVhHNWdZR0JjYm1CZ1lDSjkgLS0+XG5cbmBgYHJcbmBgYHJcbmdncGxvdChkYXRhPWRmX0xBVVIsIGFlcyh4PW1vbnRoLCB5PWhlaWdodF9tbSwgZ3JvdXA9IG1vbnRoKSkgK1xuICBnZW9tX2JveHBsb3QoKSsgICNzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXM9YyhcXGZvcmVzdGdyZWVuXFwsXFxvcmFuZ2VcXCwgXFxwdXJwbGVcXCkpK1xuICB0aGVtZV9idygpICsgIHRoZW1lKHBhbmVsLmdyaWQubWFqb3IgPSBlbGVtZW50X2JsYW5rKCksIHBhbmVsLmdyaWQubWlub3IgPSBlbGVtZW50X2JsYW5rKCkpKyBcbiAgI3RoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gOTAsIHZqdXN0ID0gMC41LCBoanVzdD0xKSkrXG4gIGxhYnModGl0bGU9XFwgTGF1cmVsIEhvbGxvdyBNZWFuIFNoZWxsIEhlaWdodCAtIDIwMjNcXCwgeCA9XFxtb250aFxcLCB5ID0gXFxNZWFuIFNoZWxsIEhlaWdodCAobW0pXFwpXG5cblxuYGBgXG5gYGBcblxuPCEtLSBybmItc291cmNlLWVuZCAtLT5cbiJ9 -->
<!-- rnb-source-begin eyJkYXRhIjoiYGBgclxuYGBgclxuZ2dwbG90KGRhdGE9ZGZfTEFVUiwgYWVzKHg9bW9udGgsIHk9aGVpZ2h0X21tLCBncm91cD0gbW9udGgpKSArXG4gIGdlb21fYm94cGxvdCgpKyAgI3NjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcz1jKFxcZm9yZXN0Z3JlZW5cXCxcXG9yYW5nZVxcLCBcXHB1cnBsZVxcKSkrXG4gIHRoZW1lX2J3KCkgKyAgdGhlbWUocGFuZWwuZ3JpZC5tYWpvciA9IGVsZW1lbnRfYmxhbmsoKSwgcGFuZWwuZ3JpZC5taW5vciA9IGVsZW1lbnRfYmxhbmsoKSkrIFxuICAjdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA5MCwgdmp1c3QgPSAwLjUsIGhqdXN0PTEpKStcbiAgbGFicyh0aXRsZT1cXCBMYXVyZWwgSG9sbG93IE1lYW4gU2hlbGwgSGVpZ2h0IC0gMjAyM1xcLCB4ID1cXG1vbnRoXFwsIHkgPSBcXE1lYW4gU2hlbGwgSGVpZ2h0IChtbSlcXClcblxuXG5gYGBcbmBgYCJ9 -->
```r
```r
ggplot(data=df_LAUR, aes(x=month, y=height_mm, group= month)) +
geom_boxplot()+ #scale_fill_manual(values=c(\forestgreen\,\orange\, \purple\))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title=\ Laurel Hollow Mean Shell Height - 2023\, x =\month\, y = \Mean Shell Height (mm)\)
```
```
<!-- rnb-source-end -->
<!-- rnb-output-end -->
<!-- rnb-chunk-end -->
<!-- rnb-chunk-begin -->
<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVJME52Ym1ScGRHbHZiaUJrYVhOMGNtbGlkWFJwYjI0Z1lXTnliM056SUdGc2JDQnpZVzF3YkdVZ2JXOXVkR2h6WEc1blozQnNiM1FvWkdGMFlUMWtabDlNUVZWU0xDQmhaWE1vZUQwZ1kyOXVaR2wwYVc5dUtTa2dLMXh1SUNCblpXOXRYMkpoY2lncEt5QWdJM05qWVd4bFgyWnBiR3hmYldGdWRXRnNLSFpoYkhWbGN6MWpLRndpWm05eVpYTjBaM0psWlc1Y0lpeGNJbTl5WVc1blpWd2lMQ0JjSW5CMWNuQnNaVndpS1NrclhHNGdJSFJvWlcxbFgySjNLQ2tnS3lBZ2RHaGxiV1VvY0dGdVpXd3VaM0pwWkM1dFlXcHZjaUE5SUdWc1pXMWxiblJmWW14aGJtc29LU3dnY0dGdVpXd3VaM0pwWkM1dGFXNXZjaUE5SUdWc1pXMWxiblJmWW14aGJtc29LU2tySUZ4dUlDQWpkR2hsYldVb1lYaHBjeTUwWlhoMExuZ2dQU0JsYkdWdFpXNTBYM1JsZUhRb1lXNW5iR1VnUFNBNU1Dd2dkbXAxYzNRZ1BTQXdMalVzSUdocWRYTjBQVEVwS1N0Y2JpQWdiR0ZpY3loMGFYUnNaVDFjSWt4aGRYSmxiQ0JJYjJ4c2IzY2dRbTlrZVNCRGIyNWthWFJwYjI0Z2FXNWtaWGdnTFNBeU1ESXpYQ0lzSUhnZ1BWd2lRMjl1WkdsMGFXOXVJR05oZEdWbmIzSnBlbUYwYVc5dVhDSXBLMXh1SUNBZ0lHWmhZMlYwWDNkeVlYQW9maUJ0YjI1MGFDd2djMk5oYkdWeklEMGdYQ0ptY21WbFhDSXBYRzVjYmx4dUkwMWxZVzRnUW05a2VTQmpiMjVrYVhScGIyNGdjR1Z5SUcxdmJuUm9YRzVuWjNCc2IzUW9aR0YwWVQxa1psOU1RVlZTTENCaFpYTW9lRDBnYlc5dWRHZ3NJSGs5SUdOdmJtUnBkR2x2Ymw5elkyOXlaU3dnWjNKdmRYQWdQU0J0YjI1MGFDa3BJQ3RjYmlBZ1oyVnZiVjlpYjNod2JHOTBLQ2tySUNBamMyTmhiR1ZmWm1sc2JGOXRZVzUxWVd3b2RtRnNkV1Z6UFdNb1hDSm1iM0psYzNSbmNtVmxibHdpTEZ3aWIzSmhibWRsWENJc0lGd2ljSFZ5Y0d4bFhDSXBLU3RjYmlBZ2RHaGxiV1ZmWW5jb0tTQXJJQ0IwYUdWdFpTaHdZVzVsYkM1bmNtbGtMbTFoYW05eUlEMGdaV3hsYldWdWRGOWliR0Z1YXlncExDQndZVzVsYkM1bmNtbGtMbTFwYm05eUlEMGdaV3hsYldWdWRGOWliR0Z1YXlncEtTc2dYRzRnSUd4aFluTW9kR2wwYkdVOVhDSk1ZWFZ5Wld3Z1NHOXNiRzkzSUUxbFlXNGdRbTlrZVNCRGIyNWthWFJwYjI0Z1UyTnZjbVVnTFNBeU1ESXpYQ0lzSUhnZ1BWd2liVzl1ZEdoY0lpd2dlVDBnWENJZ1kyOXVaR2wwYVc5dUlITmpiM0psSUNneExUa3BYQ0lwS3lCelkyRnNaVjk1WDNKbGRtVnljMlVvS1Z4dUlDQWpkR2hsYldVb1lYaHBjeTUwWlhoMExuZ2dQU0JsYkdWdFpXNTBYM1JsZUhRb1lXNW5iR1VnUFNBNU1Dd2dkbXAxYzNRZ1BTQXdMalVzSUdocWRYTjBQVEVwS1N0Y2JtQmdZQ0o5IC0tPlxuXG5gYGByXG4jQ29uZGl0aW9uIGRpc3RyaWJ1dGlvbiBhY3Jvc3MgYWxsIHNhbXBsZSBtb250aHNcbmdncGxvdChkYXRhPWRmX0xBVVIsIGFlcyh4PSBjb25kaXRpb24pKSArXG4gIGdlb21fYmFyKCkrICAjc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzPWMoXFxmb3Jlc3RncmVlblxcLFxcb3JhbmdlXFwsIFxccHVycGxlXFwpKStcbiAgdGhlbWVfYncoKSArICB0aGVtZShwYW5lbC5ncmlkLm1ham9yID0gZWxlbWVudF9ibGFuaygpLCBwYW5lbC5ncmlkLm1pbm9yID0gZWxlbWVudF9ibGFuaygpKSsgXG4gICN0aGVtZShheGlzLnRleHQueCA9IGVsZW1lbnRfdGV4dChhbmdsZSA9IDkwLCB2anVzdCA9IDAuNSwgaGp1c3Q9MSkpK1xuICBsYWJzKHRpdGxlPVxcTGF1cmVsIEhvbGxvdyBCb2R5IENvbmRpdGlvbiBpbmRleCAtIDIwMjNcXCwgeCA9XFxDb25kaXRpb24gY2F0ZWdvcml6YXRpb25cXCkrXG4gICAgZmFjZXRfd3JhcCh+IG1vbnRoLCBzY2FsZXMgPSBcXGZyZWVcXClcblxuXG4jTWVhbiBCb2R5IGNvbmRpdGlvbiBwZXIgbW9udGhcbmdncGxvdChkYXRhPWRmX0xBVVIsIGFlcyh4PSBtb250aCwgeT0gY29uZGl0aW9uX3Njb3JlLCBncm91cCA9IG1vbnRoKSkgK1xuICBnZW9tX2JveHBsb3QoKSsgICNzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXM9YyhcXGZvcmVzdGdyZWVuXFwsXFxvcmFuZ2VcXCwgXFxwdXJwbGVcXCkpK1xuICB0aGVtZV9idygpICsgIHRoZW1lKHBhbmVsLmdyaWQubWFqb3IgPSBlbGVtZW50X2JsYW5rKCksIHBhbmVsLmdyaWQubWlub3IgPSBlbGVtZW50X2JsYW5rKCkpKyBcbiAgbGFicyh0aXRsZT1cXExhdXJlbCBIb2xsb3cgTWVhbiBCb2R5IENvbmRpdGlvbiBTY29yZSAtIDIwMjNcXCwgeCA9XFxtb250aFxcLCB5PSBcXCBjb25kaXRpb24gc2NvcmUgKDEtOSlcXCkrIHNjYWxlX3lfcmV2ZXJzZSgpXG4gICN0aGVtZShheGlzLnRleHQueCA9IGVsZW1lbnRfdGV4dChhbmdsZSA9IDkwLCB2anVzdCA9IDAuNSwgaGp1c3Q9MSkpK1xuYGBgXG5cbjwhLS0gcm5iLXNvdXJjZS1lbmQgLS0+XG4ifQ== -->
#Condition distribution across all sample months
ggplot(data=df_LAUR, aes(x= condition)) +
geom_bar()+ #scale_fill_manual(values=c(\forestgreen\,\orange\, \purple\))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title=\Laurel Hollow Body Condition index - 2023\, x =\Condition categorization\)+
facet_wrap(~ month, scales = \free\)
#Mean Body condition per month
ggplot(data=df_LAUR, aes(x= month, y= condition_score, group = month)) +
geom_boxplot()+ #scale_fill_manual(values=c(\forestgreen\,\orange\, \purple\))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
labs(title=\Laurel Hollow Mean Body Condition Score - 2023\, x =\month\, y= \ condition score (1-9)\)+ scale_y_reverse()
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
<!-- rnb-output-end -->
<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVZR0JnY2x4dUkwTnZibVJwZEdsdmJpQmthWE4wY21saWRYUnBiMjRnWVdOeWIzTnpJR0ZzYkNCellXMXdiR1VnYlc5dWRHaHpYRzVuWjNCc2IzUW9aR0YwWVQxa1psOU1RVlZTTENCaFpYTW9lRDBnWTI5dVpHbDBhVzl1S1NrZ0sxeHVJQ0JuWlc5dFgySmhjaWdwS3lBZ0kzTmpZV3hsWDJacGJHeGZiV0Z1ZFdGc0tIWmhiSFZsY3oxaktGeGNabTl5WlhOMFozSmxaVzVjWEN4Y1hHOXlZVzVuWlZ4Y0xDQmNYSEIxY25Cc1pWeGNLU2tyWEc0Z0lIUm9aVzFsWDJKM0tDa2dLeUFnZEdobGJXVW9jR0Z1Wld3dVozSnBaQzV0WVdwdmNpQTlJR1ZzWlcxbGJuUmZZbXhoYm1zb0tTd2djR0Z1Wld3dVozSnBaQzV0YVc1dmNpQTlJR1ZzWlcxbGJuUmZZbXhoYm1zb0tTa3JJRnh1SUNBamRHaGxiV1VvWVhocGN5NTBaWGgwTG5nZ1BTQmxiR1Z0Wlc1MFgzUmxlSFFvWVc1bmJHVWdQU0E1TUN3Z2RtcDFjM1FnUFNBd0xqVXNJR2hxZFhOMFBURXBLU3RjYmlBZ2JHRmljeWgwYVhSc1pUMWNYRXhoZFhKbGJDQkliMnhzYjNjZ1FtOWtlU0JEYjI1a2FYUnBiMjRnYVc1a1pYZ2dMU0F5TURJelhGd3NJSGdnUFZ4Y1EyOXVaR2wwYVc5dUlHTmhkR1ZuYjNKcGVtRjBhVzl1WEZ3cEsxeHVJQ0FnSUdaaFkyVjBYM2R5WVhBb2ZpQnRiMjUwYUN3Z2MyTmhiR1Z6SUQwZ1hGeG1jbVZsWEZ3cFhHNWNibHh1STAxbFlXNGdRbTlrZVNCamIyNWthWFJwYjI0Z2NHVnlJRzF2Ym5Sb1hHNW5aM0JzYjNRb1pHRjBZVDFrWmw5TVFWVlNMQ0JoWlhNb2VEMGdiVzl1ZEdnc0lIazlJR052Ym1ScGRHbHZibDl6WTI5eVpTd2daM0p2ZFhBZ1BTQnRiMjUwYUNrcElDdGNiaUFnWjJWdmJWOWliM2h3Ykc5MEtDa3JJQ0FqYzJOaGJHVmZabWxzYkY5dFlXNTFZV3dvZG1Gc2RXVnpQV01vWEZ4bWIzSmxjM1JuY21WbGJseGNMRnhjYjNKaGJtZGxYRndzSUZ4Y2NIVnljR3hsWEZ3cEtTdGNiaUFnZEdobGJXVmZZbmNvS1NBcklDQjBhR1Z0WlNod1lXNWxiQzVuY21sa0xtMWhhbTl5SUQwZ1pXeGxiV1Z1ZEY5aWJHRnVheWdwTENCd1lXNWxiQzVuY21sa0xtMXBibTl5SUQwZ1pXeGxiV1Z1ZEY5aWJHRnVheWdwS1NzZ1hHNGdJR3hoWW5Nb2RHbDBiR1U5WEZ4TVlYVnlaV3dnU0c5c2JHOTNJRTFsWVc0Z1FtOWtlU0JEYjI1a2FYUnBiMjRnVTJOdmNtVWdMU0F5TURJelhGd3NJSGdnUFZ4Y2JXOXVkR2hjWEN3Z2VUMGdYRndnWTI5dVpHbDBhVzl1SUhOamIzSmxJQ2d4TFRrcFhGd3BLeUJ6WTJGc1pWOTVYM0psZG1WeWMyVW9LVnh1SUNBamRHaGxiV1VvWVhocGN5NTBaWGgwTG5nZ1BTQmxiR1Z0Wlc1MFgzUmxlSFFvWVc1bmJHVWdQU0E1TUN3Z2RtcDFjM1FnUFNBd0xqVXNJR2hxZFhOMFBURXBLU3RjYm1CZ1lGeHVZR0JnSW4wPSAtLT5cblxuYGBgclxuYGBgclxuI0NvbmRpdGlvbiBkaXN0cmlidXRpb24gYWNyb3NzIGFsbCBzYW1wbGUgbW9udGhzXG5nZ3Bsb3QoZGF0YT1kZl9MQVVSLCBhZXMoeD0gY29uZGl0aW9uKSkgK1xuICBnZW9tX2JhcigpKyAgI3NjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcz1jKFxcZm9yZXN0Z3JlZW5cXCxcXG9yYW5nZVxcLCBcXHB1cnBsZVxcKSkrXG4gIHRoZW1lX2J3KCkgKyAgdGhlbWUocGFuZWwuZ3JpZC5tYWpvciA9IGVsZW1lbnRfYmxhbmsoKSwgcGFuZWwuZ3JpZC5taW5vciA9IGVsZW1lbnRfYmxhbmsoKSkrIFxuICAjdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA5MCwgdmp1c3QgPSAwLjUsIGhqdXN0PTEpKStcbiAgbGFicyh0aXRsZT1cXExhdXJlbCBIb2xsb3cgQm9keSBDb25kaXRpb24gaW5kZXggLSAyMDIzXFwsIHggPVxcQ29uZGl0aW9uIGNhdGVnb3JpemF0aW9uXFwpK1xuICAgIGZhY2V0X3dyYXAofiBtb250aCwgc2NhbGVzID0gXFxmcmVlXFwpXG5cblxuI01lYW4gQm9keSBjb25kaXRpb24gcGVyIG1vbnRoXG5nZ3Bsb3QoZGF0YT1kZl9MQVVSLCBhZXMoeD0gbW9udGgsIHk9IGNvbmRpdGlvbl9zY29yZSwgZ3JvdXAgPSBtb250aCkpICtcbiAgZ2VvbV9ib3hwbG90KCkrICAjc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzPWMoXFxmb3Jlc3RncmVlblxcLFxcb3JhbmdlXFwsIFxccHVycGxlXFwpKStcbiAgdGhlbWVfYncoKSArICB0aGVtZShwYW5lbC5ncmlkLm1ham9yID0gZWxlbWVudF9ibGFuaygpLCBwYW5lbC5ncmlkLm1pbm9yID0gZWxlbWVudF9ibGFuaygpKSsgXG4gIGxhYnModGl0bGU9XFxMYXVyZWwgSG9sbG93IE1lYW4gQm9keSBDb25kaXRpb24gU2NvcmUgLSAyMDIzXFwsIHggPVxcbW9udGhcXCwgeT0gXFwgY29uZGl0aW9uIHNjb3JlICgxLTkpXFwpKyBzY2FsZV95X3JldmVyc2UoKVxuICAjdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA5MCwgdmp1c3QgPSAwLjUsIGhqdXN0PTEpKStcbmBgYFxuYGBgXG5cbjwhLS0gcm5iLXNvdXJjZS1lbmQgLS0+XG4ifQ== -->
<!-- rnb-source-begin eyJkYXRhIjoiYGBgclxuYGBgclxuI0NvbmRpdGlvbiBkaXN0cmlidXRpb24gYWNyb3NzIGFsbCBzYW1wbGUgbW9udGhzXG5nZ3Bsb3QoZGF0YT1kZl9MQVVSLCBhZXMoeD0gY29uZGl0aW9uKSkgK1xuICBnZW9tX2JhcigpKyAgI3NjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcz1jKFxcZm9yZXN0Z3JlZW5cXCxcXG9yYW5nZVxcLCBcXHB1cnBsZVxcKSkrXG4gIHRoZW1lX2J3KCkgKyAgdGhlbWUocGFuZWwuZ3JpZC5tYWpvciA9IGVsZW1lbnRfYmxhbmsoKSwgcGFuZWwuZ3JpZC5taW5vciA9IGVsZW1lbnRfYmxhbmsoKSkrIFxuICAjdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA5MCwgdmp1c3QgPSAwLjUsIGhqdXN0PTEpKStcbiAgbGFicyh0aXRsZT1cXExhdXJlbCBIb2xsb3cgQm9keSBDb25kaXRpb24gaW5kZXggLSAyMDIzXFwsIHggPVxcQ29uZGl0aW9uIGNhdGVnb3JpemF0aW9uXFwpK1xuICAgIGZhY2V0X3dyYXAofiBtb250aCwgc2NhbGVzID0gXFxmcmVlXFwpXG5cblxuI01lYW4gQm9keSBjb25kaXRpb24gcGVyIG1vbnRoXG5nZ3Bsb3QoZGF0YT1kZl9MQVVSLCBhZXMoeD0gbW9udGgsIHk9IGNvbmRpdGlvbl9zY29yZSwgZ3JvdXAgPSBtb250aCkpICtcbiAgZ2VvbV9ib3hwbG90KCkrICAjc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzPWMoXFxmb3Jlc3RncmVlblxcLFxcb3JhbmdlXFwsIFxccHVycGxlXFwpKStcbiAgdGhlbWVfYncoKSArICB0aGVtZShwYW5lbC5ncmlkLm1ham9yID0gZWxlbWVudF9ibGFuaygpLCBwYW5lbC5ncmlkLm1pbm9yID0gZWxlbWVudF9ibGFuaygpKSsgXG4gIGxhYnModGl0bGU9XFxMYXVyZWwgSG9sbG93IE1lYW4gQm9keSBDb25kaXRpb24gU2NvcmUgLSAyMDIzXFwsIHggPVxcbW9udGhcXCwgeT0gXFwgY29uZGl0aW9uIHNjb3JlICgxLTkpXFwpKyBzY2FsZV95X3JldmVyc2UoKVxuICAjdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA5MCwgdmp1c3QgPSAwLjUsIGhqdXN0PTEpKStcbmBgYFxuYGBgIn0= -->
```r
```r
#Condition distribution across all sample months
ggplot(data=df_LAUR, aes(x= condition)) +
geom_bar()+ #scale_fill_manual(values=c(\forestgreen\,\orange\, \purple\))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title=\Laurel Hollow Body Condition index - 2023\, x =\Condition categorization\)+
facet_wrap(~ month, scales = \free\)
#Mean Body condition per month
ggplot(data=df_LAUR, aes(x= month, y= condition_score, group = month)) +
geom_boxplot()+ #scale_fill_manual(values=c(\forestgreen\,\orange\, \purple\))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
labs(title=\Laurel Hollow Mean Body Condition Score - 2023\, x =\month\, y= \ condition score (1-9)\)+ scale_y_reverse()
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
```
```
<!-- rnb-source-end -->
<!-- rnb-output-end -->
<!-- rnb-chunk-end -->
<!-- rnb-text-begin -->
#LAUR Proportions graph Body condition scores - 2023 (only sampled one year)
<!-- rnb-text-end -->
<!-- rnb-chunk-begin -->
<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVYRzVNUVZWU1gwSkRVMTl3Y205d2IzSjBhVzl1Y3p3dElHUm1YMHhCVlZJZ0pUNGxYRzRnSUcxMWRHRjBaU2hqYjI1a2FYUnBiMjVmYzJOdmNtVmZiblZ0WlhKcFl5QTlJR0Z6TG01MWJXVnlhV01vWTI5dVpHbDBhVzl1WDNOamIzSmxLU3hqYjI1a2FYUnBiMjVmYzJOdmNtVmZZbWx1SUQwZ1kyRnpaVjkzYUdWdUtHTnZibVJwZEdsdmJsOXpZMjl5WlNBOVBTQXhJSDRnWENJeFhDSXNJR052Ym1ScGRHbHZibDl6WTI5eVpTQTlQU0F5SUg1Y0lqSmNJaXhqYjI1a2FYUnBiMjVmYzJOdmNtVWdQVDBnTXlCK0lGd2lNMXdpTENCamIyNWthWFJwYjI1ZmMyTnZjbVVnUFQwZ05DQitYQ0kwWENJc1kyOXVaR2wwYVc5dVgzTmpiM0psSUQwOUlEVWdmaUJjSWpWY0lpd2dZMjl1WkdsMGFXOXVYM05qYjNKbElEMDlJRFlnZmx3aU5sd2lMR052Ym1ScGRHbHZibDl6WTI5eVpTQTlQU0EzSUg0Z1hDSTNYQ0lzSUdOdmJtUnBkR2x2Ymw5elkyOXlaU0E5UFNBNElINWNJamhjSWl3Z1kyOXVaR2wwYVc5dVgzTmpiM0psSUQwOU9TQitYQ0k1WENJc0lGUlNWVVVnZmlCaGN5NWphR0Z5WVdOMFpYSW9ZMjl1WkdsMGFXOXVYM05qYjNKbEtTa3BJQ1UrSlZ4dUlDQm5jbTkxY0Y5aWVTaHphWFJsTENCdGIyNTBhQ3dnWTI5dVpHbDBhVzl1WDNOamIzSmxYMkpwYml3Z2VXVmhjaWtnSlQ0bElHUndiSGx5T2pwemRXMXRZWEpwYzJVb1EyOTFiblE5SUc0b0tTa2dKVDRsWEc0Z0lIVnVaM0p2ZFhBb0tTQWxQaVZjYmlBZ2JYVjBZWFJsS0ZCeWIzQnZjblJwYjI0Z1BTQkRiM1Z1ZEM5emRXMG9RMjkxYm5RcEtWeHVYRzVNUVZWU1gwSkRVMTl3Y205d2IzSjBhVzl1Y3p3dElHNWhMbTl0YVhRb1RFRlZVbDlDUTFOZmNISnZjRzl5ZEdsdmJuTXBYRzVjYmt4QlZWSmZRa05UWDNCeWIzQnZjblJwYjI1elhHNWNia0pEVTE5d2NtOXdiM0owYVc5dVgweEJWVkk4TFNCblozQnNiM1FvWkdGMFlUMU1RVlZTWDBKRFUxOXdjbTl3YjNKMGFXOXVjeXdnWVdWektIZzliVzl1ZEdnc0lIazlJRkJ5YjNCdmNuUnBiMjRzSUdacGJHdzlZMjl1WkdsMGFXOXVYM05qYjNKbFgySnBiaWtwSUN0Y2JpQWdaMlZ2YlY5aVlYSW9kMmxrZEdnZ1BTQXVOU3dnYzNSaGREMWNJbWxrWlc1MGFYUjVYQ0lzSUhCdmMybDBhVzl1SUQwZ1hDSm1hV3hzWENJcEt5QWdYRzRnSUNCMGFHVnRaVjlpZHlncElDc2dJSFJvWlcxbEtIQmhibVZzTG1keWFXUXViV0ZxYjNJZ1BTQmxiR1Z0Wlc1MFgySnNZVzVyS0Nrc0lIQmhibVZzTG1keWFXUXViV2x1YjNJZ1BTQmxiR1Z0Wlc1MFgySnNZVzVyS0NrcEt5QmNiaUFnZEdobGJXVW9ZWGhwY3k1MFpYaDBMbmdnUFNCbGJHVnRaVzUwWDNSbGVIUW9ZVzVuYkdVZ1BTQTBOU3dnZG1wMWMzUWdQU0F4TENCb2FuVnpkRDB4S1NrclhHNGdJR3hoWW5Nb2RHbDBiR1U5WENKTVlYVnlaV3dnU0c5c2JHOTNJRkJ5YjNCdmNuUnBiMjRnYjJZZ1ltOWtlU0JqYjI1a2FYUnBiMjRnYzJOdmNtVnpJQzBnTWpBeU0xd2lMQ0I0SUQxY0ltMXZiblJvWENJc0lIa2dQU0JjSWxCeWIzQnZjblJwYjI0Z2IyWWdZbTlrZVNCamIyNWthWFJwYjI0Z2MyTnZjbVZ6SUZ3aUtTc2dkR2hsYldVb1lYaHBjeTUwYVhSc1pTNTVJRDBnWld4bGJXVnVkRjkwWlhoMEtITnBlbVVnUFNCeVpXd29NUzR6S1N3Z1lXNW5iR1VnUFRrd0tTd2dZWGhwY3k1MGFYUnNaUzU0SUQwZ1pXeGxiV1Z1ZEY5MFpYaDBLSE5wZW1VZ1BTQnlaV3dvTVM0ektTd2dZVzVuYkdVZ1BTQXdLU2tyWEc0Z0lIUm9aVzFsS0dGNGFYTXVkR1Y0ZEQxbGJHVnRaVzUwWDNSbGVIUW9jMmw2WlQweE1pa3BYRzRnSUNObVlXTmxkRjkzY21Gd0tINGdlV1ZoY2lsY2JpQWpjMk5oYkdWZlptbHNiRjlpY21WM1pYSW9LU0FyWEc0Z0lDTm1ZV05sZEY5M2NtRndLSDRnYzJsMFpTbGNibHh1UWtOVFgzQnliM0J2Y25ScGIyNWZURUZWVWx4dUlDQmNibkJrWmlod1lYTjBaVEFvY0dGMGFDQTlJRndpVEdGaVgwUmhkR0ZmVkdsemMzVmxVSEp2WTJWemMybHVaeTl2ZFhSd2RYUmNJaUFzWENJdlRFRlZVbDlDUTFOZmNISnZjRzl5ZEdsdmJpNXdaR1pjSWlrc0lHaGxhV2RvZENBOUlEY3NJSGRwWkhSb0lEMGdNVE1wWEc1d2NtbHVkQ2hDUTFOZmNISnZjRzl5ZEdsdmJsOU1RVlZTS1Z4dVpHVjJMbTltWmlncFhHNWdZR0FpZlE9PSAtLT5cblxuYGBgclxuXG5MQVVSX0JDU19wcm9wb3J0aW9uczwtIGRmX0xBVVIgJT4lXG4gIG11dGF0ZShjb25kaXRpb25fc2NvcmVfbnVtZXJpYyA9IGFzLm51bWVyaWMoY29uZGl0aW9uX3Njb3JlKSxjb25kaXRpb25fc2NvcmVfYmluID0gY2FzZV93aGVuKGNvbmRpdGlvbl9zY29yZSA9PSAxIH4gXFwxXFwsIGNvbmRpdGlvbl9zY29yZSA9PSAyIH5cXDJcXCxjb25kaXRpb25fc2NvcmUgPT0gMyB+IFxcM1xcLCBjb25kaXRpb25fc2NvcmUgPT0gNCB+XFw0XFwsY29uZGl0aW9uX3Njb3JlID09IDUgfiBcXDVcXCwgY29uZGl0aW9uX3Njb3JlID09IDYgflxcNlxcLGNvbmRpdGlvbl9zY29yZSA9PSA3IH4gXFw3XFwsIGNvbmRpdGlvbl9zY29yZSA9PSA4IH5cXDhcXCwgY29uZGl0aW9uX3Njb3JlID09OSB+XFw5XFwsIFRSVUUgfiBhcy5jaGFyYWN0ZXIoY29uZGl0aW9uX3Njb3JlKSkpICU+JVxuICBncm91cF9ieShzaXRlLCBtb250aCwgY29uZGl0aW9uX3Njb3JlX2JpbiwgeWVhcikgJT4lIGRwbHlyOjpzdW1tYXJpc2UoQ291bnQ9IG4oKSkgJT4lXG4gIHVuZ3JvdXAoKSAlPiVcbiAgbXV0YXRlKFByb3BvcnRpb24gPSBDb3VudC9zdW0oQ291bnQpKVxuXG5MQVVSX0JDU19wcm9wb3J0aW9uczwtIG5hLm9taXQoTEFVUl9CQ1NfcHJvcG9ydGlvbnMpXG5cbkxBVVJfQkNTX3Byb3BvcnRpb25zXG5cbkJDU19wcm9wb3J0aW9uX0xBVVI8LSBnZ3Bsb3QoZGF0YT1MQVVSX0JDU19wcm9wb3J0aW9ucywgYWVzKHg9bW9udGgsIHk9IFByb3BvcnRpb24sIGZpbGw9Y29uZGl0aW9uX3Njb3JlX2JpbikpICtcbiAgZ2VvbV9iYXIod2lkdGggPSAuNSwgc3RhdD1cXGlkZW50aXR5XFwsIHBvc2l0aW9uID0gXFxmaWxsXFwpKyAgXG4gICB0aGVtZV9idygpICsgIHRoZW1lKHBhbmVsLmdyaWQubWFqb3IgPSBlbGVtZW50X2JsYW5rKCksIHBhbmVsLmdyaWQubWlub3IgPSBlbGVtZW50X2JsYW5rKCkpKyBcbiAgdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA0NSwgdmp1c3QgPSAxLCBoanVzdD0xKSkrXG4gIGxhYnModGl0bGU9XFxMYXVyZWwgSG9sbG93IFByb3BvcnRpb24gb2YgYm9keSBjb25kaXRpb24gc2NvcmVzIC0gMjAyM1xcLCB4ID1cXG1vbnRoXFwsIHkgPSBcXFByb3BvcnRpb24gb2YgYm9keSBjb25kaXRpb24gc2NvcmVzIFxcKSsgdGhlbWUoYXhpcy50aXRsZS55ID0gZWxlbWVudF90ZXh0KHNpemUgPSByZWwoMS4zKSwgYW5nbGUgPTkwKSwgYXhpcy50aXRsZS54ID0gZWxlbWVudF90ZXh0KHNpemUgPSByZWwoMS4zKSwgYW5nbGUgPSAwKSkrXG4gIHRoZW1lKGF4aXMudGV4dD1lbGVtZW50X3RleHQoc2l6ZT0xMikpXG4gICNmYWNldF93cmFwKH4geWVhcilcbiAjc2NhbGVfZmlsbF9icmV3ZXIoKSArXG4gICNmYWNldF93cmFwKH4gc2l0ZSlcblxuQkNTX3Byb3BvcnRpb25fTEFVUlxuICBcbnBkZihwYXN0ZTAocGF0aCA9IFxcTGFiX0RhdGFfVGlzc3VlUHJvY2Vzc2luZy9vdXRwdXRcXCAsXFwvTEFVUl9CQ1NfcHJvcG9ydGlvbi5wZGZcXCksIGhlaWdodCA9IDcsIHdpZHRoID0gMTMpXG5wcmludChCQ1NfcHJvcG9ydGlvbl9MQVVSKVxuZGV2Lm9mZigpXG5gYGBcblxuPCEtLSBybmItc291cmNlLWVuZCAtLT5cbiJ9 -->
LAUR_BCS_proportions<- df_LAUR %>%
mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ \1\, condition_score == 2 ~\2\,condition_score == 3 ~ \3\, condition_score == 4 ~\4\,condition_score == 5 ~ \5\, condition_score == 6 ~\6\,condition_score == 7 ~ \7\, condition_score == 8 ~\8\, condition_score ==9 ~\9\, TRUE ~ as.character(condition_score))) %>%
group_by(site, month, condition_score_bin, year) %>% dplyr::summarise(Count= n()) %>%
ungroup() %>%
mutate(Proportion = Count/sum(Count))
LAUR_BCS_proportions<- na.omit(LAUR_BCS_proportions)
LAUR_BCS_proportions
BCS_proportion_LAUR<- ggplot(data=LAUR_BCS_proportions, aes(x=month, y= Proportion, fill=condition_score_bin)) +
geom_bar(width = .5, stat=\identity\, position = \fill\)+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
labs(title=\Laurel Hollow Proportion of body condition scores - 2023\, x =\month\, y = \Proportion of body condition scores \)+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
theme(axis.text=element_text(size=12))
#facet_wrap(~ year)
#scale_fill_brewer() +
#facet_wrap(~ site)
BCS_proportion_LAUR
pdf(paste0(path = \Lab_Data_TissueProcessing/output\ ,\/LAUR_BCS_proportion.pdf\), height = 7, width = 13)
print(BCS_proportion_LAUR)
dev.off()
<!-- rnb-output-end -->
<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVZR0JnY2x4dVhHNU1RVlZTWDBKRFUxOXdjbTl3YjNKMGFXOXVjend0SUdSbVgweEJWVklnSlQ0bFhHNGdJRzExZEdGMFpTaGpiMjVrYVhScGIyNWZjMk52Y21WZmJuVnRaWEpwWXlBOUlHRnpMbTUxYldWeWFXTW9ZMjl1WkdsMGFXOXVYM05qYjNKbEtTeGpiMjVrYVhScGIyNWZjMk52Y21WZlltbHVJRDBnWTJGelpWOTNhR1Z1S0dOdmJtUnBkR2x2Ymw5elkyOXlaU0E5UFNBeElINGdYRnd4WEZ3c0lHTnZibVJwZEdsdmJsOXpZMjl5WlNBOVBTQXlJSDVjWERKY1hDeGpiMjVrYVhScGIyNWZjMk52Y21VZ1BUMGdNeUIrSUZ4Y00xeGNMQ0JqYjI1a2FYUnBiMjVmYzJOdmNtVWdQVDBnTkNCK1hGdzBYRndzWTI5dVpHbDBhVzl1WDNOamIzSmxJRDA5SURVZ2ZpQmNYRFZjWEN3Z1kyOXVaR2wwYVc5dVgzTmpiM0psSUQwOUlEWWdmbHhjTmx4Y0xHTnZibVJwZEdsdmJsOXpZMjl5WlNBOVBTQTNJSDRnWEZ3M1hGd3NJR052Ym1ScGRHbHZibDl6WTI5eVpTQTlQU0E0SUg1Y1hEaGNYQ3dnWTI5dVpHbDBhVzl1WDNOamIzSmxJRDA5T1NCK1hGdzVYRndzSUZSU1ZVVWdmaUJoY3k1amFHRnlZV04wWlhJb1kyOXVaR2wwYVc5dVgzTmpiM0psS1NrcElDVStKVnh1SUNCbmNtOTFjRjlpZVNoemFYUmxMQ0J0YjI1MGFDd2dZMjl1WkdsMGFXOXVYM05qYjNKbFgySnBiaXdnZVdWaGNpa2dKVDRsSUdSd2JIbHlPanB6ZFcxdFlYSnBjMlVvUTI5MWJuUTlJRzRvS1NrZ0pUNGxYRzRnSUhWdVozSnZkWEFvS1NBbFBpVmNiaUFnYlhWMFlYUmxLRkJ5YjNCdmNuUnBiMjRnUFNCRGIzVnVkQzl6ZFcwb1EyOTFiblFwS1Z4dVhHNU1RVlZTWDBKRFUxOXdjbTl3YjNKMGFXOXVjend0SUc1aExtOXRhWFFvVEVGVlVsOUNRMU5mY0hKdmNHOXlkR2x2Ym5NcFhHNWNia3hCVlZKZlFrTlRYM0J5YjNCdmNuUnBiMjV6WEc1Y2JrSkRVMTl3Y205d2IzSjBhVzl1WDB4QlZWSThMU0JuWjNCc2IzUW9aR0YwWVQxTVFWVlNYMEpEVTE5d2NtOXdiM0owYVc5dWN5d2dZV1Z6S0hnOWJXOXVkR2dzSUhrOUlGQnliM0J2Y25ScGIyNHNJR1pwYkd3OVkyOXVaR2wwYVc5dVgzTmpiM0psWDJKcGJpa3BJQ3RjYmlBZ1oyVnZiVjlpWVhJb2QybGtkR2dnUFNBdU5Td2djM1JoZEQxY1hHbGtaVzUwYVhSNVhGd3NJSEJ2YzJsMGFXOXVJRDBnWEZ4bWFXeHNYRndwS3lBZ1hHNGdJQ0IwYUdWdFpWOWlkeWdwSUNzZ0lIUm9aVzFsS0hCaGJtVnNMbWR5YVdRdWJXRnFiM0lnUFNCbGJHVnRaVzUwWDJKc1lXNXJLQ2tzSUhCaGJtVnNMbWR5YVdRdWJXbHViM0lnUFNCbGJHVnRaVzUwWDJKc1lXNXJLQ2twS3lCY2JpQWdkR2hsYldVb1lYaHBjeTUwWlhoMExuZ2dQU0JsYkdWdFpXNTBYM1JsZUhRb1lXNW5iR1VnUFNBME5Td2dkbXAxYzNRZ1BTQXhMQ0JvYW5WemREMHhLU2tyWEc0Z0lHeGhZbk1vZEdsMGJHVTlYRnhNWVhWeVpXd2dTRzlzYkc5M0lGQnliM0J2Y25ScGIyNGdiMllnWW05a2VTQmpiMjVrYVhScGIyNGdjMk52Y21WeklDMGdNakF5TTF4Y0xDQjRJRDFjWEcxdmJuUm9YRndzSUhrZ1BTQmNYRkJ5YjNCdmNuUnBiMjRnYjJZZ1ltOWtlU0JqYjI1a2FYUnBiMjRnYzJOdmNtVnpJRnhjS1NzZ2RHaGxiV1VvWVhocGN5NTBhWFJzWlM1NUlEMGdaV3hsYldWdWRGOTBaWGgwS0hOcGVtVWdQU0J5Wld3b01TNHpLU3dnWVc1bmJHVWdQVGt3S1N3Z1lYaHBjeTUwYVhSc1pTNTRJRDBnWld4bGJXVnVkRjkwWlhoMEtITnBlbVVnUFNCeVpXd29NUzR6S1N3Z1lXNW5iR1VnUFNBd0tTa3JYRzRnSUhSb1pXMWxLR0Y0YVhNdWRHVjRkRDFsYkdWdFpXNTBYM1JsZUhRb2MybDZaVDB4TWlrcFhHNGdJQ05tWVdObGRGOTNjbUZ3S0g0Z2VXVmhjaWxjYmlBamMyTmhiR1ZmWm1sc2JGOWljbVYzWlhJb0tTQXJYRzRnSUNObVlXTmxkRjkzY21Gd0tINGdjMmwwWlNsY2JseHVRa05UWDNCeWIzQnZjblJwYjI1ZlRFRlZVbHh1SUNCY2JuQmtaaWh3WVhOMFpUQW9jR0YwYUNBOUlGeGNUR0ZpWDBSaGRHRmZWR2x6YzNWbFVISnZZMlZ6YzJsdVp5OXZkWFJ3ZFhSY1hDQXNYRnd2VEVGVlVsOUNRMU5mY0hKdmNHOXlkR2x2Ymk1d1pHWmNYQ2tzSUdobGFXZG9kQ0E5SURjc0lIZHBaSFJvSUQwZ01UTXBYRzV3Y21sdWRDaENRMU5mY0hKdmNHOXlkR2x2Ymw5TVFWVlNLVnh1WkdWMkxtOW1aaWdwWEc1Z1lHQmNibUJnWUNKOSAtLT5cblxuYGBgclxuYGBgclxuXG5MQVVSX0JDU19wcm9wb3J0aW9uczwtIGRmX0xBVVIgJT4lXG4gIG11dGF0ZShjb25kaXRpb25fc2NvcmVfbnVtZXJpYyA9IGFzLm51bWVyaWMoY29uZGl0aW9uX3Njb3JlKSxjb25kaXRpb25fc2NvcmVfYmluID0gY2FzZV93aGVuKGNvbmRpdGlvbl9zY29yZSA9PSAxIH4gXFwxXFwsIGNvbmRpdGlvbl9zY29yZSA9PSAyIH5cXDJcXCxjb25kaXRpb25fc2NvcmUgPT0gMyB+IFxcM1xcLCBjb25kaXRpb25fc2NvcmUgPT0gNCB+XFw0XFwsY29uZGl0aW9uX3Njb3JlID09IDUgfiBcXDVcXCwgY29uZGl0aW9uX3Njb3JlID09IDYgflxcNlxcLGNvbmRpdGlvbl9zY29yZSA9PSA3IH4gXFw3XFwsIGNvbmRpdGlvbl9zY29yZSA9PSA4IH5cXDhcXCwgY29uZGl0aW9uX3Njb3JlID09OSB+XFw5XFwsIFRSVUUgfiBhcy5jaGFyYWN0ZXIoY29uZGl0aW9uX3Njb3JlKSkpICU+JVxuICBncm91cF9ieShzaXRlLCBtb250aCwgY29uZGl0aW9uX3Njb3JlX2JpbiwgeWVhcikgJT4lIGRwbHlyOjpzdW1tYXJpc2UoQ291bnQ9IG4oKSkgJT4lXG4gIHVuZ3JvdXAoKSAlPiVcbiAgbXV0YXRlKFByb3BvcnRpb24gPSBDb3VudC9zdW0oQ291bnQpKVxuXG5MQVVSX0JDU19wcm9wb3J0aW9uczwtIG5hLm9taXQoTEFVUl9CQ1NfcHJvcG9ydGlvbnMpXG5cbkxBVVJfQkNTX3Byb3BvcnRpb25zXG5cbkJDU19wcm9wb3J0aW9uX0xBVVI8LSBnZ3Bsb3QoZGF0YT1MQVVSX0JDU19wcm9wb3J0aW9ucywgYWVzKHg9bW9udGgsIHk9IFByb3BvcnRpb24sIGZpbGw9Y29uZGl0aW9uX3Njb3JlX2JpbikpICtcbiAgZ2VvbV9iYXIod2lkdGggPSAuNSwgc3RhdD1cXGlkZW50aXR5XFwsIHBvc2l0aW9uID0gXFxmaWxsXFwpKyAgXG4gICB0aGVtZV9idygpICsgIHRoZW1lKHBhbmVsLmdyaWQubWFqb3IgPSBlbGVtZW50X2JsYW5rKCksIHBhbmVsLmdyaWQubWlub3IgPSBlbGVtZW50X2JsYW5rKCkpKyBcbiAgdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA0NSwgdmp1c3QgPSAxLCBoanVzdD0xKSkrXG4gIGxhYnModGl0bGU9XFxMYXVyZWwgSG9sbG93IFByb3BvcnRpb24gb2YgYm9keSBjb25kaXRpb24gc2NvcmVzIC0gMjAyM1xcLCB4ID1cXG1vbnRoXFwsIHkgPSBcXFByb3BvcnRpb24gb2YgYm9keSBjb25kaXRpb24gc2NvcmVzIFxcKSsgdGhlbWUoYXhpcy50aXRsZS55ID0gZWxlbWVudF90ZXh0KHNpemUgPSByZWwoMS4zKSwgYW5nbGUgPTkwKSwgYXhpcy50aXRsZS54ID0gZWxlbWVudF90ZXh0KHNpemUgPSByZWwoMS4zKSwgYW5nbGUgPSAwKSkrXG4gIHRoZW1lKGF4aXMudGV4dD1lbGVtZW50X3RleHQoc2l6ZT0xMikpXG4gICNmYWNldF93cmFwKH4geWVhcilcbiAjc2NhbGVfZmlsbF9icmV3ZXIoKSArXG4gICNmYWNldF93cmFwKH4gc2l0ZSlcblxuQkNTX3Byb3BvcnRpb25fTEFVUlxuICBcbnBkZihwYXN0ZTAocGF0aCA9IFxcTGFiX0RhdGFfVGlzc3VlUHJvY2Vzc2luZy9vdXRwdXRcXCAsXFwvTEFVUl9CQ1NfcHJvcG9ydGlvbi5wZGZcXCksIGhlaWdodCA9IDcsIHdpZHRoID0gMTMpXG5wcmludChCQ1NfcHJvcG9ydGlvbl9MQVVSKVxuZGV2Lm9mZigpXG5gYGBcbmBgYFxuXG48IS0tIHJuYi1zb3VyY2UtZW5kIC0tPlxuIn0= -->
<!-- rnb-source-begin eyJkYXRhIjoiYGBgclxuYGBgclxuXG5MQVVSX0JDU19wcm9wb3J0aW9uczwtIGRmX0xBVVIgJT4lXG4gIG11dGF0ZShjb25kaXRpb25fc2NvcmVfbnVtZXJpYyA9IGFzLm51bWVyaWMoY29uZGl0aW9uX3Njb3JlKSxjb25kaXRpb25fc2NvcmVfYmluID0gY2FzZV93aGVuKGNvbmRpdGlvbl9zY29yZSA9PSAxIH4gXFwxXFwsIGNvbmRpdGlvbl9zY29yZSA9PSAyIH5cXDJcXCxjb25kaXRpb25fc2NvcmUgPT0gMyB+IFxcM1xcLCBjb25kaXRpb25fc2NvcmUgPT0gNCB+XFw0XFwsY29uZGl0aW9uX3Njb3JlID09IDUgfiBcXDVcXCwgY29uZGl0aW9uX3Njb3JlID09IDYgflxcNlxcLGNvbmRpdGlvbl9zY29yZSA9PSA3IH4gXFw3XFwsIGNvbmRpdGlvbl9zY29yZSA9PSA4IH5cXDhcXCwgY29uZGl0aW9uX3Njb3JlID09OSB+XFw5XFwsIFRSVUUgfiBhcy5jaGFyYWN0ZXIoY29uZGl0aW9uX3Njb3JlKSkpICU+JVxuICBncm91cF9ieShzaXRlLCBtb250aCwgY29uZGl0aW9uX3Njb3JlX2JpbiwgeWVhcikgJT4lIGRwbHlyOjpzdW1tYXJpc2UoQ291bnQ9IG4oKSkgJT4lXG4gIHVuZ3JvdXAoKSAlPiVcbiAgbXV0YXRlKFByb3BvcnRpb24gPSBDb3VudC9zdW0oQ291bnQpKVxuXG5MQVVSX0JDU19wcm9wb3J0aW9uczwtIG5hLm9taXQoTEFVUl9CQ1NfcHJvcG9ydGlvbnMpXG5cbkxBVVJfQkNTX3Byb3BvcnRpb25zXG5cbkJDU19wcm9wb3J0aW9uX0xBVVI8LSBnZ3Bsb3QoZGF0YT1MQVVSX0JDU19wcm9wb3J0aW9ucywgYWVzKHg9bW9udGgsIHk9IFByb3BvcnRpb24sIGZpbGw9Y29uZGl0aW9uX3Njb3JlX2JpbikpICtcbiAgZ2VvbV9iYXIod2lkdGggPSAuNSwgc3RhdD1cXGlkZW50aXR5XFwsIHBvc2l0aW9uID0gXFxmaWxsXFwpKyAgXG4gICB0aGVtZV9idygpICsgIHRoZW1lKHBhbmVsLmdyaWQubWFqb3IgPSBlbGVtZW50X2JsYW5rKCksIHBhbmVsLmdyaWQubWlub3IgPSBlbGVtZW50X2JsYW5rKCkpKyBcbiAgdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA0NSwgdmp1c3QgPSAxLCBoanVzdD0xKSkrXG4gIGxhYnModGl0bGU9XFxMYXVyZWwgSG9sbG93IFByb3BvcnRpb24gb2YgYm9keSBjb25kaXRpb24gc2NvcmVzIC0gMjAyM1xcLCB4ID1cXG1vbnRoXFwsIHkgPSBcXFByb3BvcnRpb24gb2YgYm9keSBjb25kaXRpb24gc2NvcmVzIFxcKSsgdGhlbWUoYXhpcy50aXRsZS55ID0gZWxlbWVudF90ZXh0KHNpemUgPSByZWwoMS4zKSwgYW5nbGUgPTkwKSwgYXhpcy50aXRsZS54ID0gZWxlbWVudF90ZXh0KHNpemUgPSByZWwoMS4zKSwgYW5nbGUgPSAwKSkrXG4gIHRoZW1lKGF4aXMudGV4dD1lbGVtZW50X3RleHQoc2l6ZT0xMikpXG4gICNmYWNldF93cmFwKH4geWVhcilcbiAjc2NhbGVfZmlsbF9icmV3ZXIoKSArXG4gICNmYWNldF93cmFwKH4gc2l0ZSlcblxuQkNTX3Byb3BvcnRpb25fTEFVUlxuICBcbnBkZihwYXN0ZTAocGF0aCA9IFxcTGFiX0RhdGFfVGlzc3VlUHJvY2Vzc2luZy9vdXRwdXRcXCAsXFwvTEFVUl9CQ1NfcHJvcG9ydGlvbi5wZGZcXCksIGhlaWdodCA9IDcsIHdpZHRoID0gMTMpXG5wcmludChCQ1NfcHJvcG9ydGlvbl9MQVVSKVxuZGV2Lm9mZigpXG5gYGBcbmBgYCJ9 -->
```r
```r
LAUR_BCS_proportions<- df_LAUR %>%
mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ \1\, condition_score == 2 ~\2\,condition_score == 3 ~ \3\, condition_score == 4 ~\4\,condition_score == 5 ~ \5\, condition_score == 6 ~\6\,condition_score == 7 ~ \7\, condition_score == 8 ~\8\, condition_score ==9 ~\9\, TRUE ~ as.character(condition_score))) %>%
group_by(site, month, condition_score_bin, year) %>% dplyr::summarise(Count= n()) %>%
ungroup() %>%
mutate(Proportion = Count/sum(Count))
LAUR_BCS_proportions<- na.omit(LAUR_BCS_proportions)
LAUR_BCS_proportions
BCS_proportion_LAUR<- ggplot(data=LAUR_BCS_proportions, aes(x=month, y= Proportion, fill=condition_score_bin)) +
geom_bar(width = .5, stat=\identity\, position = \fill\)+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
labs(title=\Laurel Hollow Proportion of body condition scores - 2023\, x =\month\, y = \Proportion of body condition scores \)+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
theme(axis.text=element_text(size=12))
#facet_wrap(~ year)
#scale_fill_brewer() +
#facet_wrap(~ site)
BCS_proportion_LAUR
pdf(paste0(path = \Lab_Data_TissueProcessing/output\ ,\/LAUR_BCS_proportion.pdf\), height = 7, width = 13)
print(BCS_proportion_LAUR)
dev.off()
```
```
<!-- rnb-source-end -->
<!-- rnb-output-end -->
<!-- rnb-chunk-end -->
<!-- rnb-text-begin -->
# % of scores >3 at LAUR june - november
<!-- rnb-text-end -->
<!-- rnb-chunk-begin -->
<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVYRzRqY21WdGIzWnBibWNnTURneU0wWkZUa05mTWpnZ1ltTnpJSGRoY3lCT1FWeHVaR1pmVEVGVlVpNGdQQzBnWkdaZlRFRlZVaVUrSldacGJIUmxjaWdoY205M1gyNTFiV0psY2lncElDVnBiaVVnWXlneUtTbGNibHh1VEVGVlVsOXdjbVZqWlc1MFgyZHlaV0YwWlhKZk15QThMU0JrWmw5TVFWVlNMaUFsUGlWY2JpQWdaSEJzZVhJNk9tZHliM1Z3WDJKNUtHMXZiblJvTENCemFYUmxMQ0I1WldGeUtTQWxQaVZjYmlBZ1pIQnNlWEk2T25OMWJXMWhjbWx6WlNoUVpYSmpaVzUwWVdkbElEMGdiV1ZoYmloamIyNWthWFJwYjI1ZmMyTnZjbVVnUEQwZ015a3FNVEF3S1Z4dVRFRlZVbDl3Y21WalpXNTBYMmR5WldGMFpYSmZNMXh1WEc1TVFWVlNYMEpEVTE5d1pYSmpaVzUwWVdkbElEd3RJRXhCVlZKZmNISmxZMlZ1ZEY5bmNtVmhkR1Z5WHpNZ0pUNGxYRzRnSUdkbmNHeHZkQ2hoWlhNb2VDQTlJRzF2Ym5Sb0xDQjVJRDBnVUdWeVkyVnVkR0ZuWlNrcElDdGNiaUFnSTJkbGIyMWZZbUZ5S0hkcFpIUm9JRDBnTUM0MUxDQnpkR0YwSUQwZ1hDSnBaR1Z1ZEdsMGVWd2lMQ0J3YjNOcGRHbHZiaUE5SUZ3aVptbHNiRndpS1NBclhHNGdJQ0FnWjJWdmJWOWpiMndvS1N0Y2JpQWdkR2hsYldWZlluY29LU0FyWEc0Z0lIUm9aVzFsS0Z4dUlDQWdJSEJoYm1Wc0xtZHlhV1F1YldGcWIzSWdQU0JsYkdWdFpXNTBYMkpzWVc1cktDa3NYRzRnSUNBZ2NHRnVaV3d1WjNKcFpDNXRhVzV2Y2lBOUlHVnNaVzFsYm5SZllteGhibXNvS1N4Y2JpQWdJQ0JoZUdsekxuUmxlSFF1ZUNBOUlHVnNaVzFsYm5SZmRHVjRkQ2hoYm1kc1pTQTlJRFExTENCMmFuVnpkQ0E5SURFc0lHaHFkWE4wSUQwZ01Ta3NYRzRnSUNBZ1lYaHBjeTUwYVhSc1pTNTVJRDBnWld4bGJXVnVkRjkwWlhoMEtITnBlbVVnUFNCeVpXd29NUzR6S1N3Z1lXNW5iR1VnUFNBNU1Da3NYRzRnSUNBZ1lYaHBjeTUwYVhSc1pTNTRJRDBnWld4bGJXVnVkRjkwWlhoMEtITnBlbVVnUFNCeVpXd29NUzR6S1N3Z1lXNW5iR1VnUFNBd0tTeGNiaUFnSUNCaGVHbHpMblJsZUhRZ1BTQmxiR1Z0Wlc1MFgzUmxlSFFvYzJsNlpTQTlJREV5S1Z4dUlDQXBJQ3RjYmlBZ2JHRmljeWhjYmlBZ0lDQjBhWFJzWlNBOUlGd2lSMjlzWkNCVGRHRnlJRUpsWVdOb0lDVWdiMllnWW05a2VTQmpiMjVrYVhScGIyNGdjMk52Y21WeklENDlJRE5jSWl4Y2JpQWdJQ0I0SUQwZ1hDSnRiMjUwYUZ3aUxGeHVJQ0FnSUhrZ1BTQmNJbEJsY21ObGJuUmhaMlVnYjJZZ1ltOWtlU0JqYjI1a2FYUnBiMjRnYzJOdmNtVnpJRDQ5SUROY0lseHVJQ0FwSUNNclptRmpaWFJmZDNKaGNDaCtlV1ZoY2lsY2JpQWdJeXNnYzJOaGJHVmZlRjlqYjI1MGFXNTFiM1Z6S0NCaWNtVmhhM01nUFNCelpYRW9OU3d4TWl3Z1lua2dQVEVwSUNrclhHNGdJQ041YkdsdEtEQXNOakFwS3lCY2JpQWdJM05qWVd4bFgzbGZZMjl1ZEdsdWRXOTFjeWhzYVcxcGRITWdQU0JqS0RBc01UQXdLU3dnWW5KbFlXdHpJRDBnYzJWeEtEQXNNVEF3TENCaWVTQTlJREV3S1NrZ1hHNWNiaUFnSUNCY2JreEJWVkpmUWtOVFgzQmxjbU5sYm5SaFoyVmNibHh1Y0dSbUtIQmhjM1JsTUNod1lYUm9JRDBnWENKTVlXSmZSR0YwWVY5VWFYTnpkV1ZRY205alpYTnphVzVuTDI5MWRIQjFkRndpTENCY0lpOU1RVlZTWDBKRFUxOXdaWEpqWlc1MFlXZGxMbkJrWmx3aUtTeG9aV2xuYUhRZ1BTQTNMQ0IzYVdSMGFDQTlJREV6S1Z4dWNISnBiblFvVEVGVlVsOUNRMU5mY0dWeVkyVnVkR0ZuWlNsY2JtUmxkaTV2Wm1Zb0tTQmNibHh1WEc1Y2JtQmdZQ0o5IC0tPlxuXG5gYGByXG5cbiNyZW1vdmluZyAwODIzRkVOQ18yOCBiY3Mgd2FzIE5BXG5kZl9MQVVSLiA8LSBkZl9MQVVSJT4lZmlsdGVyKCFyb3dfbnVtYmVyKCkgJWluJSBjKDIpKVxuXG5MQVVSX3ByZWNlbnRfZ3JlYXRlcl8zIDwtIGRmX0xBVVIuICU+JVxuICBkcGx5cjo6Z3JvdXBfYnkobW9udGgsIHNpdGUsIHllYXIpICU+JVxuICBkcGx5cjo6c3VtbWFyaXNlKFBlcmNlbnRhZ2UgPSBtZWFuKGNvbmRpdGlvbl9zY29yZSA8PSAzKSoxMDApXG5MQVVSX3ByZWNlbnRfZ3JlYXRlcl8zXG5cbkxBVVJfQkNTX3BlcmNlbnRhZ2UgPC0gTEFVUl9wcmVjZW50X2dyZWF0ZXJfMyAlPiVcbiAgZ2dwbG90KGFlcyh4ID0gbW9udGgsIHkgPSBQZXJjZW50YWdlKSkgK1xuICAjZ2VvbV9iYXIod2lkdGggPSAwLjUsIHN0YXQgPSBcXGlkZW50aXR5XFwsIHBvc2l0aW9uID0gXFxmaWxsXFwpICtcbiAgICBnZW9tX2NvbCgpK1xuICB0aGVtZV9idygpICtcbiAgdGhlbWUoXG4gICAgcGFuZWwuZ3JpZC5tYWpvciA9IGVsZW1lbnRfYmxhbmsoKSxcbiAgICBwYW5lbC5ncmlkLm1pbm9yID0gZWxlbWVudF9ibGFuaygpLFxuICAgIGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gNDUsIHZqdXN0ID0gMSwgaGp1c3QgPSAxKSxcbiAgICBheGlzLnRpdGxlLnkgPSBlbGVtZW50X3RleHQoc2l6ZSA9IHJlbCgxLjMpLCBhbmdsZSA9IDkwKSxcbiAgICBheGlzLnRpdGxlLnggPSBlbGVtZW50X3RleHQoc2l6ZSA9IHJlbCgxLjMpLCBhbmdsZSA9IDApLFxuICAgIGF4aXMudGV4dCA9IGVsZW1lbnRfdGV4dChzaXplID0gMTIpXG4gICkgK1xuICBsYWJzKFxuICAgIHRpdGxlID0gXFxHb2xkIFN0YXIgQmVhY2ggJSBvZiBib2R5IGNvbmRpdGlvbiBzY29yZXMgPj0gM1xcLFxuICAgIHggPSBcXG1vbnRoXFwsXG4gICAgeSA9IFxcUGVyY2VudGFnZSBvZiBib2R5IGNvbmRpdGlvbiBzY29yZXMgPj0gM1xcXG4gICkgIytmYWNldF93cmFwKH55ZWFyKVxuICAjKyBzY2FsZV94X2NvbnRpbnVvdXMoIGJyZWFrcyA9IHNlcSg1LDEyLCBieSA9MSkgKStcbiAgI3lsaW0oMCw2MCkrIFxuICAjc2NhbGVfeV9jb250aW51b3VzKGxpbWl0cyA9IGMoMCwxMDApLCBicmVha3MgPSBzZXEoMCwxMDAsIGJ5ID0gMTApKSBcblxuICAgIFxuTEFVUl9CQ1NfcGVyY2VudGFnZVxuXG5wZGYocGFzdGUwKHBhdGggPSBcXExhYl9EYXRhX1Rpc3N1ZVByb2Nlc3Npbmcvb3V0cHV0XFwsIFxcL0xBVVJfQkNTX3BlcmNlbnRhZ2UucGRmXFwpLGhlaWdodCA9IDcsIHdpZHRoID0gMTMpXG5wcmludChMQVVSX0JDU19wZXJjZW50YWdlKVxuZGV2Lm9mZigpIFxuXG5cblxuYGBgXG5cbjwhLS0gcm5iLXNvdXJjZS1lbmQgLS0+XG4ifQ== -->
#removing 0823FENC_28 bcs was NA
df_LAUR. <- df_LAUR%>%filter(!row_number() %in% c(2))
LAUR_precent_greater_3 <- df_LAUR. %>%
dplyr::group_by(month, site, year) %>%
dplyr::summarise(Percentage = mean(condition_score <= 3)*100)
LAUR_precent_greater_3
LAUR_BCS_percentage <- LAUR_precent_greater_3 %>%
ggplot(aes(x = month, y = Percentage)) +
#geom_bar(width = 0.5, stat = \identity\, position = \fill\) +
geom_col()+
theme_bw() +
theme(
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1),
axis.title.y = element_text(size = rel(1.3), angle = 90),
axis.title.x = element_text(size = rel(1.3), angle = 0),
axis.text = element_text(size = 12)
) +
labs(
title = \Gold Star Beach % of body condition scores >= 3\,
x = \month\,
y = \Percentage of body condition scores >= 3\
) #+facet_wrap(~year)
#+ scale_x_continuous( breaks = seq(5,12, by =1) )+
#ylim(0,60)+
#scale_y_continuous(limits = c(0,100), breaks = seq(0,100, by = 10))
LAUR_BCS_percentage
pdf(paste0(path = \Lab_Data_TissueProcessing/output\, \/LAUR_BCS_percentage.pdf\),height = 7, width = 13)
print(LAUR_BCS_percentage)
dev.off()
<!-- rnb-output-end -->
<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVZR0JnY2x4dVhHNGpjbVZ0YjNacGJtY2dNRGd5TTBaRlRrTmZNamdnWW1OeklIZGhjeUJPUVZ4dVpHWmZURUZWVWk0Z1BDMGdaR1pmVEVGVlVpVStKV1pwYkhSbGNpZ2hjbTkzWDI1MWJXSmxjaWdwSUNWcGJpVWdZeWd5S1NsY2JseHVURUZWVWw5d2NtVmpaVzUwWDJkeVpXRjBaWEpmTXlBOExTQmtabDlNUVZWU0xpQWxQaVZjYmlBZ1pIQnNlWEk2T21keWIzVndYMko1S0cxdmJuUm9MQ0J6YVhSbExDQjVaV0Z5S1NBbFBpVmNiaUFnWkhCc2VYSTZPbk4xYlcxaGNtbHpaU2hRWlhKalpXNTBZV2RsSUQwZ2JXVmhiaWhqYjI1a2FYUnBiMjVmYzJOdmNtVWdQRDBnTXlrcU1UQXdLVnh1VEVGVlVsOXdjbVZqWlc1MFgyZHlaV0YwWlhKZk0xeHVYRzVNUVZWU1gwSkRVMTl3WlhKalpXNTBZV2RsSUR3dElFeEJWVkpmY0hKbFkyVnVkRjluY21WaGRHVnlYek1nSlQ0bFhHNGdJR2RuY0d4dmRDaGhaWE1vZUNBOUlHMXZiblJvTENCNUlEMGdVR1Z5WTJWdWRHRm5aU2twSUN0Y2JpQWdJMmRsYjIxZlltRnlLSGRwWkhSb0lEMGdNQzQxTENCemRHRjBJRDBnWEZ4cFpHVnVkR2wwZVZ4Y0xDQndiM05wZEdsdmJpQTlJRnhjWm1sc2JGeGNLU0FyWEc0Z0lDQWdaMlZ2YlY5amIyd29LU3RjYmlBZ2RHaGxiV1ZmWW5jb0tTQXJYRzRnSUhSb1pXMWxLRnh1SUNBZ0lIQmhibVZzTG1keWFXUXViV0ZxYjNJZ1BTQmxiR1Z0Wlc1MFgySnNZVzVyS0Nrc1hHNGdJQ0FnY0dGdVpXd3VaM0pwWkM1dGFXNXZjaUE5SUdWc1pXMWxiblJmWW14aGJtc29LU3hjYmlBZ0lDQmhlR2x6TG5SbGVIUXVlQ0E5SUdWc1pXMWxiblJmZEdWNGRDaGhibWRzWlNBOUlEUTFMQ0IyYW5WemRDQTlJREVzSUdocWRYTjBJRDBnTVNrc1hHNGdJQ0FnWVhocGN5NTBhWFJzWlM1NUlEMGdaV3hsYldWdWRGOTBaWGgwS0hOcGVtVWdQU0J5Wld3b01TNHpLU3dnWVc1bmJHVWdQU0E1TUNrc1hHNGdJQ0FnWVhocGN5NTBhWFJzWlM1NElEMGdaV3hsYldWdWRGOTBaWGgwS0hOcGVtVWdQU0J5Wld3b01TNHpLU3dnWVc1bmJHVWdQU0F3S1N4Y2JpQWdJQ0JoZUdsekxuUmxlSFFnUFNCbGJHVnRaVzUwWDNSbGVIUW9jMmw2WlNBOUlERXlLVnh1SUNBcElDdGNiaUFnYkdGaWN5aGNiaUFnSUNCMGFYUnNaU0E5SUZ4Y1IyOXNaQ0JUZEdGeUlFSmxZV05vSUNVZ2IyWWdZbTlrZVNCamIyNWthWFJwYjI0Z2MyTnZjbVZ6SUQ0OUlETmNYQ3hjYmlBZ0lDQjRJRDBnWEZ4dGIyNTBhRnhjTEZ4dUlDQWdJSGtnUFNCY1hGQmxjbU5sYm5SaFoyVWdiMllnWW05a2VTQmpiMjVrYVhScGIyNGdjMk52Y21WeklENDlJRE5jWEZ4dUlDQXBJQ01yWm1GalpYUmZkM0poY0NoK2VXVmhjaWxjYmlBZ0l5c2djMk5oYkdWZmVGOWpiMjUwYVc1MWIzVnpLQ0JpY21WaGEzTWdQU0J6WlhFb05Td3hNaXdnWW5rZ1BURXBJQ2tyWEc0Z0lDTjViR2x0S0RBc05qQXBLeUJjYmlBZ0kzTmpZV3hsWDNsZlkyOXVkR2x1ZFc5MWN5aHNhVzFwZEhNZ1BTQmpLREFzTVRBd0tTd2dZbkpsWVd0eklEMGdjMlZ4S0RBc01UQXdMQ0JpZVNBOUlERXdLU2tnWEc1Y2JpQWdJQ0JjYmt4QlZWSmZRa05UWDNCbGNtTmxiblJoWjJWY2JseHVjR1JtS0hCaGMzUmxNQ2h3WVhSb0lEMGdYRnhNWVdKZlJHRjBZVjlVYVhOemRXVlFjbTlqWlhOemFXNW5MMjkxZEhCMWRGeGNMQ0JjWEM5TVFWVlNYMEpEVTE5d1pYSmpaVzUwWVdkbExuQmtabHhjS1N4b1pXbG5hSFFnUFNBM0xDQjNhV1IwYUNBOUlERXpLVnh1Y0hKcGJuUW9URUZWVWw5Q1ExTmZjR1Z5WTJWdWRHRm5aU2xjYm1SbGRpNXZabVlvS1NCY2JseHVYRzVjYm1CZ1lGeHVZR0JnSW4wPSAtLT5cblxuYGBgclxuYGBgclxuXG4jcmVtb3ZpbmcgMDgyM0ZFTkNfMjggYmNzIHdhcyBOQVxuZGZfTEFVUi4gPC0gZGZfTEFVUiU+JWZpbHRlcighcm93X251bWJlcigpICVpbiUgYygyKSlcblxuTEFVUl9wcmVjZW50X2dyZWF0ZXJfMyA8LSBkZl9MQVVSLiAlPiVcbiAgZHBseXI6Omdyb3VwX2J5KG1vbnRoLCBzaXRlLCB5ZWFyKSAlPiVcbiAgZHBseXI6OnN1bW1hcmlzZShQZXJjZW50YWdlID0gbWVhbihjb25kaXRpb25fc2NvcmUgPD0gMykqMTAwKVxuTEFVUl9wcmVjZW50X2dyZWF0ZXJfM1xuXG5MQVVSX0JDU19wZXJjZW50YWdlIDwtIExBVVJfcHJlY2VudF9ncmVhdGVyXzMgJT4lXG4gIGdncGxvdChhZXMoeCA9IG1vbnRoLCB5ID0gUGVyY2VudGFnZSkpICtcbiAgI2dlb21fYmFyKHdpZHRoID0gMC41LCBzdGF0ID0gXFxpZGVudGl0eVxcLCBwb3NpdGlvbiA9IFxcZmlsbFxcKSArXG4gICAgZ2VvbV9jb2woKStcbiAgdGhlbWVfYncoKSArXG4gIHRoZW1lKFxuICAgIHBhbmVsLmdyaWQubWFqb3IgPSBlbGVtZW50X2JsYW5rKCksXG4gICAgcGFuZWwuZ3JpZC5taW5vciA9IGVsZW1lbnRfYmxhbmsoKSxcbiAgICBheGlzLnRleHQueCA9IGVsZW1lbnRfdGV4dChhbmdsZSA9IDQ1LCB2anVzdCA9IDEsIGhqdXN0ID0gMSksXG4gICAgYXhpcy50aXRsZS55ID0gZWxlbWVudF90ZXh0KHNpemUgPSByZWwoMS4zKSwgYW5nbGUgPSA5MCksXG4gICAgYXhpcy50aXRsZS54ID0gZWxlbWVudF90ZXh0KHNpemUgPSByZWwoMS4zKSwgYW5nbGUgPSAwKSxcbiAgICBheGlzLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDEyKVxuICApICtcbiAgbGFicyhcbiAgICB0aXRsZSA9IFxcR29sZCBTdGFyIEJlYWNoICUgb2YgYm9keSBjb25kaXRpb24gc2NvcmVzID49IDNcXCxcbiAgICB4ID0gXFxtb250aFxcLFxuICAgIHkgPSBcXFBlcmNlbnRhZ2Ugb2YgYm9keSBjb25kaXRpb24gc2NvcmVzID49IDNcXFxuICApICMrZmFjZXRfd3JhcCh+eWVhcilcbiAgIysgc2NhbGVfeF9jb250aW51b3VzKCBicmVha3MgPSBzZXEoNSwxMiwgYnkgPTEpICkrXG4gICN5bGltKDAsNjApKyBcbiAgI3NjYWxlX3lfY29udGludW91cyhsaW1pdHMgPSBjKDAsMTAwKSwgYnJlYWtzID0gc2VxKDAsMTAwLCBieSA9IDEwKSkgXG5cbiAgICBcbkxBVVJfQkNTX3BlcmNlbnRhZ2VcblxucGRmKHBhc3RlMChwYXRoID0gXFxMYWJfRGF0YV9UaXNzdWVQcm9jZXNzaW5nL291dHB1dFxcLCBcXC9MQVVSX0JDU19wZXJjZW50YWdlLnBkZlxcKSxoZWlnaHQgPSA3LCB3aWR0aCA9IDEzKVxucHJpbnQoTEFVUl9CQ1NfcGVyY2VudGFnZSlcbmRldi5vZmYoKSBcblxuXG5cbmBgYFxuYGBgXG5cbjwhLS0gcm5iLXNvdXJjZS1lbmQgLS0+XG4ifQ== -->
<!-- rnb-source-begin eyJkYXRhIjoiYGBgclxuYGBgclxuXG4jcmVtb3ZpbmcgMDgyM0ZFTkNfMjggYmNzIHdhcyBOQVxuZGZfTEFVUi4gPC0gZGZfTEFVUiU+JWZpbHRlcighcm93X251bWJlcigpICVpbiUgYygyKSlcblxuTEFVUl9wcmVjZW50X2dyZWF0ZXJfMyA8LSBkZl9MQVVSLiAlPiVcbiAgZHBseXI6Omdyb3VwX2J5KG1vbnRoLCBzaXRlLCB5ZWFyKSAlPiVcbiAgZHBseXI6OnN1bW1hcmlzZShQZXJjZW50YWdlID0gbWVhbihjb25kaXRpb25fc2NvcmUgPD0gMykqMTAwKVxuTEFVUl9wcmVjZW50X2dyZWF0ZXJfM1xuXG5MQVVSX0JDU19wZXJjZW50YWdlIDwtIExBVVJfcHJlY2VudF9ncmVhdGVyXzMgJT4lXG4gIGdncGxvdChhZXMoeCA9IG1vbnRoLCB5ID0gUGVyY2VudGFnZSkpICtcbiAgI2dlb21fYmFyKHdpZHRoID0gMC41LCBzdGF0ID0gXFxpZGVudGl0eVxcLCBwb3NpdGlvbiA9IFxcZmlsbFxcKSArXG4gICAgZ2VvbV9jb2woKStcbiAgdGhlbWVfYncoKSArXG4gIHRoZW1lKFxuICAgIHBhbmVsLmdyaWQubWFqb3IgPSBlbGVtZW50X2JsYW5rKCksXG4gICAgcGFuZWwuZ3JpZC5taW5vciA9IGVsZW1lbnRfYmxhbmsoKSxcbiAgICBheGlzLnRleHQueCA9IGVsZW1lbnRfdGV4dChhbmdsZSA9IDQ1LCB2anVzdCA9IDEsIGhqdXN0ID0gMSksXG4gICAgYXhpcy50aXRsZS55ID0gZWxlbWVudF90ZXh0KHNpemUgPSByZWwoMS4zKSwgYW5nbGUgPSA5MCksXG4gICAgYXhpcy50aXRsZS54ID0gZWxlbWVudF90ZXh0KHNpemUgPSByZWwoMS4zKSwgYW5nbGUgPSAwKSxcbiAgICBheGlzLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDEyKVxuICApICtcbiAgbGFicyhcbiAgICB0aXRsZSA9IFxcR29sZCBTdGFyIEJlYWNoICUgb2YgYm9keSBjb25kaXRpb24gc2NvcmVzID49IDNcXCxcbiAgICB4ID0gXFxtb250aFxcLFxuICAgIHkgPSBcXFBlcmNlbnRhZ2Ugb2YgYm9keSBjb25kaXRpb24gc2NvcmVzID49IDNcXFxuICApICMrZmFjZXRfd3JhcCh+eWVhcilcbiAgIysgc2NhbGVfeF9jb250aW51b3VzKCBicmVha3MgPSBzZXEoNSwxMiwgYnkgPTEpICkrXG4gICN5bGltKDAsNjApKyBcbiAgI3NjYWxlX3lfY29udGludW91cyhsaW1pdHMgPSBjKDAsMTAwKSwgYnJlYWtzID0gc2VxKDAsMTAwLCBieSA9IDEwKSkgXG5cbiAgICBcbkxBVVJfQkNTX3BlcmNlbnRhZ2VcblxucGRmKHBhc3RlMChwYXRoID0gXFxMYWJfRGF0YV9UaXNzdWVQcm9jZXNzaW5nL291dHB1dFxcLCBcXC9MQVVSX0JDU19wZXJjZW50YWdlLnBkZlxcKSxoZWlnaHQgPSA3LCB3aWR0aCA9IDEzKVxucHJpbnQoTEFVUl9CQ1NfcGVyY2VudGFnZSlcbmRldi5vZmYoKSBcblxuXG5cbmBgYFxuYGBgIn0= -->
```r
```r
#removing 0823FENC_28 bcs was NA
df_LAUR. <- df_LAUR%>%filter(!row_number() %in% c(2))
LAUR_precent_greater_3 <- df_LAUR. %>%
dplyr::group_by(month, site, year) %>%
dplyr::summarise(Percentage = mean(condition_score <= 3)*100)
LAUR_precent_greater_3
LAUR_BCS_percentage <- LAUR_precent_greater_3 %>%
ggplot(aes(x = month, y = Percentage)) +
#geom_bar(width = 0.5, stat = \identity\, position = \fill\) +
geom_col()+
theme_bw() +
theme(
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1),
axis.title.y = element_text(size = rel(1.3), angle = 90),
axis.title.x = element_text(size = rel(1.3), angle = 0),
axis.text = element_text(size = 12)
) +
labs(
title = \Gold Star Beach % of body condition scores >= 3\,
x = \month\,
y = \Percentage of body condition scores >= 3\
) #+facet_wrap(~year)
#+ scale_x_continuous( breaks = seq(5,12, by =1) )+
#ylim(0,60)+
#scale_y_continuous(limits = c(0,100), breaks = seq(0,100, by = 10))
LAUR_BCS_percentage
pdf(paste0(path = \Lab_Data_TissueProcessing/output\, \/LAUR_BCS_percentage.pdf\),height = 7, width = 13)
print(LAUR_BCS_percentage)
dev.off()
```
```
<!-- rnb-source-end -->
<!-- rnb-output-end -->
<!-- rnb-chunk-end -->
<!-- rnb-text-begin -->
<!-- rnb-text-end -->
<!-- rnb-chunk-begin -->
<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVJME52WkdsdVp5Qm5jbUZ3YUhNZ1ptOXlJRXhKVXlCRGIyNW1aWEpsYm1ObElDMGdiMjVzZVNBeU1ESXpJR1JoZEdGY2JtUmhkR0ZmTWpBeU15QThMU0JrWVhSaFgyRnNiQ1UrSlNBZ1ptbHNkR1Z5S0NGNVpXRnlQVDFjSWpJd01qUmNJaWxjYmx4dVpHWmZRa05UWDNCeWIzQnZjblJwYjI1elBDMGdaR0YwWVY4eU1ESXpJQ1UrSlZ4dUlDQnRkWFJoZEdVb1kyOXVaR2wwYVc5dVgzTmpiM0psWDI1MWJXVnlhV01nUFNCaGN5NXVkVzFsY21saktHTnZibVJwZEdsdmJsOXpZMjl5WlNrc1kyOXVaR2wwYVc5dVgzTmpiM0psWDJKcGJpQTlJR05oYzJWZmQyaGxiaWhqYjI1a2FYUnBiMjVmYzJOdmNtVWdQVDBnTVNCK0lGd2lNVndpTENCamIyNWthWFJwYjI1ZmMyTnZjbVVnUFQwZ01pQitYQ0l5WENJc1kyOXVaR2wwYVc5dVgzTmpiM0psSUQwOUlETWdmaUJjSWpOY0lpd2dZMjl1WkdsMGFXOXVYM05qYjNKbElEMDlJRFFnZmx3aU5Gd2lMR052Ym1ScGRHbHZibDl6WTI5eVpTQTlQU0ExSUg0Z1hDSTFYQ0lzSUdOdmJtUnBkR2x2Ymw5elkyOXlaU0E5UFNBMklINWNJalpjSWl4amIyNWthWFJwYjI1ZmMyTnZjbVVnUFQwZ055QitJRndpTjF3aUxDQmpiMjVrYVhScGIyNWZjMk52Y21VZ1BUMGdPQ0IrWENJNFhDSXNJR052Ym1ScGRHbHZibDl6WTI5eVpTQTlQVGtnZmx3aU9Wd2lMQ0JVVWxWRklINGdZWE11WTJoaGNtRmpkR1Z5S0dOdmJtUnBkR2x2Ymw5elkyOXlaU2twS1NBbFBpVmNiaUFnWjNKdmRYQmZZbmtvYzJsMFpTd2diVzl1ZEdnc0lHTnZibVJwZEdsdmJsOXpZMjl5WlY5aWFXNHBJQ1UrSlNCa2NHeDVjaUE2T25OMWJXMWhjbWx6WlNoRGIzVnVkRDBnYmlncEtTQWxQaVZjYmlBZ2RXNW5jbTkxY0NncElDVStKVnh1SUNCdGRYUmhkR1VvVUhKdmNHOXlkR2x2YmlBOUlFTnZkVzUwTDNOMWJTaERiM1Z1ZENrcFhHNWNibVJtWDBKRFUxOXdjbTl3YjNKMGFXOXVjend0SUc1aExtOXRhWFFvWkdaZlFrTlRYM0J5YjNCdmNuUnBiMjV6S1Z4dVhHNWtabDlDUTFOZmNISnZjRzl5ZEdsdmJuTWtiVzl1ZEdnZ1BDMGdabUZqZEc5eUtHUm1YMEpEVTE5d2NtOXdiM0owYVc5dWN5UnRiMjUwYUN3Z2JHVjJaV3h6SUQwZ1l5aGNJalZjSWl4Y0lqWmNJaXhjSWpkY0lpd2dYQ0k0WENJc0lGd2lPVndpTENCY0lqRXdYQ0lzSUZ3aU1URmNJaWtzWEc0Z0lDQWdJQ0FnSUd4aFltVnNjejFqS0Z3aVRXRjVYQ0lzSUZ3aVNuVnVaVndpTENCY0lrcDFiSGxjSWl3Z1hDSkJkV2RjSWl3Z1hDSlRaWEIwWENJc0lGd2lUMk4wWENJc0lGd2lUbTkyWENJcEtWeHVaR1pmUWtOVFgzQnliM0J2Y25ScGIyNXpKSE5wZEdVZ1BDMGdabUZqZEc5eUtHUm1YMEpEVTE5d2NtOXdiM0owYVc5dWN5UnphWFJsTENCc1pYWmxiSE1nUFNCaktGd2lRVk5JUTF3aUxGd2lSa1ZPUTF3aUxGd2lSMDlNUkZ3aUxDQmNJa3hCVlZKY0lpa3NYRzRnSUNBZ0lDQWdJR3hoWW1Wc2N6MWpLRndpUVhOb0lFTnlaV1ZyWENJc0lGd2lSbVZ1WTJVZ1EzSmxaV3RjSWl3Z1hDSkhiMnhrSUZOMFlYSmNJaXdnWENKTVlYVnlaV3dnU0c5c2JHOTNYQ0lwS1Z4dVhHNWlZM011Y0hKdmNDNHlNREl6SUR3dElHZG5jR3h2ZENoa1lYUmhQV1JtWDBKRFUxOXdjbTl3YjNKMGFXOXVjeXdnWVdWektIZzliVzl1ZEdnc0lIazlJRkJ5YjNCdmNuUnBiMjRzSUdacGJHdzlZMjl1WkdsMGFXOXVYM05qYjNKbFgySnBiaWtwSUN0Y2JpQWdaMlZ2YlY5aVlYSW9kMmxrZEdnZ1BTQXVOU3dnYzNSaGREMWNJbWxrWlc1MGFYUjVYQ0lzSUhCdmMybDBhVzl1SUQwZ1hDSm1hV3hzWENJc0lHTnZiRzkxY2lBOUlGd2lZbXhoWTJ0Y0lpa3JJQ0JjYmlBZ0lIUm9aVzFsWDJKM0tDa2dLeUFnZEdobGJXVW9jR0Z1Wld3dVozSnBaQzV0YVc1dmNpQTlJR1ZzWlcxbGJuUmZZbXhoYm1zb0tTa3JJRnh1SUNCMGFHVnRaU2hoZUdsekxuUmxlSFF1ZUNBOUlHVnNaVzFsYm5SZmRHVjRkQ2hoYm1kc1pTQTlJRFExTENCMmFuVnpkQ0E5SURFc0lHaHFkWE4wUFRFc0lITnBlbVVnUFNCeVpXd29MamtwS1NrclhHNGdJSFJvWlcxbEtIQmhibVZzTG1keWFXUXViV0ZxYjNJdWVDQTlaV3hsYldWdWRGOWliR0Z1YXlncEtTdGNiaUFnYkdGaWN5aDBhWFJzWlQxY0lrSnZaSGtnWTI5dVpHbDBhVzl1SUhOamIzSmxjeUJwYmlBeU1ESXpYQ0lzSUhnZ1BWd2lUVzl1ZEdoY0lpd2dlU0E5SUZ3aVVISnZjRzl5ZEdsdmJpQnZaaUJpYjJSNUlHTnZibVJwZEdsdmJpQnpZMjl5WlhNZ1hDSXBLeUIwYUdWdFpTaGhlR2x6TG5ScGRHeGxMbmtnUFNCbGJHVnRaVzUwWDNSbGVIUW9jMmw2WlNBOUlISmxiQ2d4TGpJcExDQmhibWRzWlNBOU9UQXBMQ0JoZUdsekxuUnBkR3hsTG5nZ1BTQmxiR1Z0Wlc1MFgzUmxlSFFvYzJsNlpTQTlJSEpsYkNneExqSXBMQ0JoYm1kc1pTQTlJREFwS1NzZ2MyTmhiR1ZmWm1sc2JGOWljbVYzWlhJb2NHRnNaWFIwWlNBOUlGd2lRbXgxWlhOY0lpa2dLMXh1SUNCMGFHVnRaU2h6ZEhKcGNDNTBaWGgwSUQwZ1pXeGxiV1Z1ZEY5MFpYaDBLSE5wZW1VZ1BTQXhOU2twSzF4dUlDQjBhR1Z0WlNBb2JHVm5aVzVrTG5ScGRHeGxJRDBnWld4bGJXVnVkRjkwWlhoMEtITnBlbVVnUFNBeE9Da3BLMXh1SUNCMGFHVnRaU2gwYVhSc1pTQTlJR1ZzWlcxbGJuUmZkR1Y0ZENoemFYcGxJRDBnTVRjcEtTdGNiaUFnZEdobGJXVW9ZWGhwY3k1MFpYaDBQV1ZzWlcxbGJuUmZkR1Y0ZENoemFYcGxQVEUzS1NrclhHNGdJR3hoWW5Nb1ptbHNiQ0E5WENKQ2IyUjVJRU52Ym1ScGRHbHZiaUJUWTI5eVpWd2lLU3RjYmlBZ1ptRmpaWFJmZDNKaGNDaCtJSE5wZEdVZ0tWeHVZbU56TG5CeWIzQXVNakF5TTF4dVhHNGpjR1JtS0hCaGMzUmxNQ2h3WVhSb0lEMGdYQ0pNWVdKZlJHRjBZVjlVYVhOemRXVlFjbTlqWlhOemFXNW5MMjkxZEhCMWRGd2lJQ3hjSWk5Q1ExTmZjSEp2Y0c5eWRHbHZibDh5TURJekxuQmtabHdpS1N3Z2FHVnBaMmgwSUQwZ055d2dkMmxrZEdnZ1BTQXhNeWxjYmlOd2NtbHVkQ2hpWTNNdWNISnZjQzR5TURJektWeHVJMlJsZGk1dlptWW9LVnh1WEc0alFrTlRYM0J5YjNCdmNuUnBiMjVmWVd4c0szTmpZV3hsWDJacGJHeGZiV0Z1ZFdGc0tIWmhiSFZsY3lBOUlHTW9YQ0lqTURBelF6TXdYQ0lzWENJak1ERTJOalZGWENJc1hDSWpPREJEUkVNeFhDSXNYQ0lqUXpkRlFVVTFYQ0lzWENJalJqWkZPRU16WENJc1hDSWpSRVpETWpkRVhDSXNYQ0lqUWtZNE1USkVYQ0lzSUZ3aUl6aEROVEV3UVZ3aUxGd2lJelUwTXpBd05Wd2lJQ2twWEc0amMyTmhiR1ZmWm1sc2JGOXRZVzUxWVd3b2RtRnNkV1Z6SUQwZ1l5aGNJaU13TURORE16QmNJaXhjSWlNd01UWTJOVVZjSWl4Y0lpTTRNRU5FUXpGY0lpeGNJaU5ETjBWQlJUVmNJaXhjSW1GNmRYSmxNVndpTEZ3aWMyeGhkR1ZuY21GNU1Wd2lMRndpYzJ4aGRHVm5jbUY1TTF3aUxDQmNJbk5zWVhSbFozSmhlVFJjSWl3Z1hDSm5jbUY1TWpWY0lpQXBLU0JjYmx4dVppQThMU0JtZFc1amRHbHZiaWh3WVd3cElHSnlaWGRsY2k1d1lXd29ZbkpsZDJWeUxuQmhiQzVwYm1adlczQmhiQ3dnWENKdFlYaGpiMnh2Y25OY0lsMHNJSEJoYkNsY2JpaGpiMnh6SUR3dElHWW9YQ0paYkVkdVFuVmNJaWtwWEc1Z1lHQWlmUT09IC0tPlxuXG5gYGByXG4jQ29kaW5nIGdyYXBocyBmb3IgTElTIENvbmZlcmVuY2UgLSBvbmx5IDIwMjMgZGF0YVxuZGF0YV8yMDIzIDwtIGRhdGFfYWxsJT4lICBmaWx0ZXIoIXllYXI9PVxcMjAyNFxcKVxuXG5kZl9CQ1NfcHJvcG9ydGlvbnM8LSBkYXRhXzIwMjMgJT4lXG4gIG11dGF0ZShjb25kaXRpb25fc2NvcmVfbnVtZXJpYyA9IGFzLm51bWVyaWMoY29uZGl0aW9uX3Njb3JlKSxjb25kaXRpb25fc2NvcmVfYmluID0gY2FzZV93aGVuKGNvbmRpdGlvbl9zY29yZSA9PSAxIH4gXFwxXFwsIGNvbmRpdGlvbl9zY29yZSA9PSAyIH5cXDJcXCxjb25kaXRpb25fc2NvcmUgPT0gMyB+IFxcM1xcLCBjb25kaXRpb25fc2NvcmUgPT0gNCB+XFw0XFwsY29uZGl0aW9uX3Njb3JlID09IDUgfiBcXDVcXCwgY29uZGl0aW9uX3Njb3JlID09IDYgflxcNlxcLGNvbmRpdGlvbl9zY29yZSA9PSA3IH4gXFw3XFwsIGNvbmRpdGlvbl9zY29yZSA9PSA4IH5cXDhcXCwgY29uZGl0aW9uX3Njb3JlID09OSB+XFw5XFwsIFRSVUUgfiBhcy5jaGFyYWN0ZXIoY29uZGl0aW9uX3Njb3JlKSkpICU+JVxuICBncm91cF9ieShzaXRlLCBtb250aCwgY29uZGl0aW9uX3Njb3JlX2JpbikgJT4lIGRwbHlyIDo6c3VtbWFyaXNlKENvdW50PSBuKCkpICU+JVxuICB1bmdyb3VwKCkgJT4lXG4gIG11dGF0ZShQcm9wb3J0aW9uID0gQ291bnQvc3VtKENvdW50KSlcblxuZGZfQkNTX3Byb3BvcnRpb25zPC0gbmEub21pdChkZl9CQ1NfcHJvcG9ydGlvbnMpXG5cbmRmX0JDU19wcm9wb3J0aW9ucyRtb250aCA8LSBmYWN0b3IoZGZfQkNTX3Byb3BvcnRpb25zJG1vbnRoLCBsZXZlbHMgPSBjKFxcNVxcLFxcNlxcLFxcN1xcLCBcXDhcXCwgXFw5XFwsIFxcMTBcXCwgXFwxMVxcKSxcbiAgICAgICAgbGFiZWxzPWMoXFxNYXlcXCwgXFxKdW5lXFwsIFxcSnVseVxcLCBcXEF1Z1xcLCBcXFNlcHRcXCwgXFxPY3RcXCwgXFxOb3ZcXCkpXG5kZl9CQ1NfcHJvcG9ydGlvbnMkc2l0ZSA8LSBmYWN0b3IoZGZfQkNTX3Byb3BvcnRpb25zJHNpdGUsIGxldmVscyA9IGMoXFxBU0hDXFwsXFxGRU5DXFwsXFxHT0xEXFwsIFxcTEFVUlxcKSxcbiAgICAgICAgbGFiZWxzPWMoXFxBc2ggQ3JlZWtcXCwgXFxGZW5jZSBDcmVla1xcLCBcXEdvbGQgU3RhclxcLCBcXExhdXJlbCBIb2xsb3dcXCkpXG5cbmJjcy5wcm9wLjIwMjMgPC0gZ2dwbG90KGRhdGE9ZGZfQkNTX3Byb3BvcnRpb25zLCBhZXMoeD1tb250aCwgeT0gUHJvcG9ydGlvbiwgZmlsbD1jb25kaXRpb25fc2NvcmVfYmluKSkgK1xuICBnZW9tX2Jhcih3aWR0aCA9IC41LCBzdGF0PVxcaWRlbnRpdHlcXCwgcG9zaXRpb24gPSBcXGZpbGxcXCwgY29sb3VyID0gXFxibGFja1xcKSsgIFxuICAgdGhlbWVfYncoKSArICB0aGVtZShwYW5lbC5ncmlkLm1pbm9yID0gZWxlbWVudF9ibGFuaygpKSsgXG4gIHRoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gNDUsIHZqdXN0ID0gMSwgaGp1c3Q9MSwgc2l6ZSA9IHJlbCguOSkpKStcbiAgdGhlbWUocGFuZWwuZ3JpZC5tYWpvci54ID1lbGVtZW50X2JsYW5rKCkpK1xuICBsYWJzKHRpdGxlPVxcQm9keSBjb25kaXRpb24gc2NvcmVzIGluIDIwMjNcXCwgeCA9XFxNb250aFxcLCB5ID0gXFxQcm9wb3J0aW9uIG9mIGJvZHkgY29uZGl0aW9uIHNjb3JlcyBcXCkrIHRoZW1lKGF4aXMudGl0bGUueSA9IGVsZW1lbnRfdGV4dChzaXplID0gcmVsKDEuMiksIGFuZ2xlID05MCksIGF4aXMudGl0bGUueCA9IGVsZW1lbnRfdGV4dChzaXplID0gcmVsKDEuMiksIGFuZ2xlID0gMCkpKyBzY2FsZV9maWxsX2JyZXdlcihwYWxldHRlID0gXFxCbHVlc1xcKSArXG4gIHRoZW1lKHN0cmlwLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSkrXG4gIHRoZW1lIChsZWdlbmQudGl0bGUgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE4KSkrXG4gIHRoZW1lKHRpdGxlID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNykpK1xuICB0aGVtZShheGlzLnRleHQ9ZWxlbWVudF90ZXh0KHNpemU9MTcpKStcbiAgbGFicyhmaWxsID1cXEJvZHkgQ29uZGl0aW9uIFNjb3JlXFwpK1xuICBmYWNldF93cmFwKH4gc2l0ZSApXG5iY3MucHJvcC4yMDIzXG5cbiNwZGYocGFzdGUwKHBhdGggPSBcXExhYl9EYXRhX1Rpc3N1ZVByb2Nlc3Npbmcvb3V0cHV0XFwgLFxcL0JDU19wcm9wb3J0aW9uXzIwMjMucGRmXFwpLCBoZWlnaHQgPSA3LCB3aWR0aCA9IDEzKVxuI3ByaW50KGJjcy5wcm9wLjIwMjMpXG4jZGV2Lm9mZigpXG5cbiNCQ1NfcHJvcG9ydGlvbl9hbGwrc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gYyhcXCMwMDNDMzBcXCxcXCMwMTY2NUVcXCxcXCM4MENEQzFcXCxcXCNDN0VBRTVcXCxcXCNGNkU4QzNcXCxcXCNERkMyN0RcXCxcXCNCRjgxMkRcXCwgXFwjOEM1MTBBXFwsXFwjNTQzMDA1XFwgKSlcbiNzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSBjKFxcIzAwM0MzMFxcLFxcIzAxNjY1RVxcLFxcIzgwQ0RDMVxcLFxcI0M3RUFFNVxcLFxcYXp1cmUxXFwsXFxzbGF0ZWdyYXkxXFwsXFxzbGF0ZWdyYXkzXFwsIFxcc2xhdGVncmF5NFxcLCBcXGdyYXkyNVxcICkpIFxuXG5mIDwtIGZ1bmN0aW9uKHBhbCkgYnJld2VyLnBhbChicmV3ZXIucGFsLmluZm9bcGFsLCBcXG1heGNvbG9yc1xcXSwgcGFsKVxuKGNvbHMgPC0gZihcXFlsR25CdVxcKSlcbmBgYFxuXG48IS0tIHJuYi1zb3VyY2UtZW5kIC0tPlxuIn0= -->
#Coding graphs for LIS Conference - only 2023 data
data_2023 <- data_all%>% filter(!year==\2024\)
df_BCS_proportions<- data_2023 %>%
mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ \1\, condition_score == 2 ~\2\,condition_score == 3 ~ \3\, condition_score == 4 ~\4\,condition_score == 5 ~ \5\, condition_score == 6 ~\6\,condition_score == 7 ~ \7\, condition_score == 8 ~\8\, condition_score ==9 ~\9\, TRUE ~ as.character(condition_score))) %>%
group_by(site, month, condition_score_bin) %>% dplyr ::summarise(Count= n()) %>%
ungroup() %>%
mutate(Proportion = Count/sum(Count))
df_BCS_proportions<- na.omit(df_BCS_proportions)
df_BCS_proportions$month <- factor(df_BCS_proportions$month, levels = c(\5\,\6\,\7\, \8\, \9\, \10\, \11\),
labels=c(\May\, \June\, \July\, \Aug\, \Sept\, \Oct\, \Nov\))
df_BCS_proportions$site <- factor(df_BCS_proportions$site, levels = c(\ASHC\,\FENC\,\GOLD\, \LAUR\),
labels=c(\Ash Creek\, \Fence Creek\, \Gold Star\, \Laurel Hollow\))
bcs.prop.2023 <- ggplot(data=df_BCS_proportions, aes(x=month, y= Proportion, fill=condition_score_bin)) +
geom_bar(width = .5, stat=\identity\, position = \fill\, colour = \black\)+
theme_bw() + theme(panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1, size = rel(.9)))+
theme(panel.grid.major.x =element_blank())+
labs(title=\Body condition scores in 2023\, x =\Month\, y = \Proportion of body condition scores \)+ theme(axis.title.y = element_text(size = rel(1.2), angle =90), axis.title.x = element_text(size = rel(1.2), angle = 0))+ scale_fill_brewer(palette = \Blues\) +
theme(strip.text = element_text(size = 15))+
theme (legend.title = element_text(size = 18))+
theme(title = element_text(size = 17))+
theme(axis.text=element_text(size=17))+
labs(fill =\Body Condition Score\)+
facet_wrap(~ site )
bcs.prop.2023
#pdf(paste0(path = \Lab_Data_TissueProcessing/output\ ,\/BCS_proportion_2023.pdf\), height = 7, width = 13)
#print(bcs.prop.2023)
#dev.off()
#BCS_proportion_all+scale_fill_manual(values = c(\#003C30\,\#01665E\,\#80CDC1\,\#C7EAE5\,\#F6E8C3\,\#DFC27D\,\#BF812D\, \#8C510A\,\#543005\ ))
#scale_fill_manual(values = c(\#003C30\,\#01665E\,\#80CDC1\,\#C7EAE5\,\azure1\,\slategray1\,\slategray3\, \slategray4\, \gray25\ ))
f <- function(pal) brewer.pal(brewer.pal.info[pal, \maxcolors\], pal)
(cols <- f(\YlGnBu\))
<!-- rnb-output-end -->
<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVZR0JnY2x4dUkwTnZaR2x1WnlCbmNtRndhSE1nWm05eUlFeEpVeUJEYjI1bVpYSmxibU5sSUMwZ2IyNXNlU0F5TURJeklHUmhkR0ZjYm1SaGRHRmZNakF5TXlBOExTQmtZWFJoWDJGc2JDVStKU0FnWm1sc2RHVnlLQ0Y1WldGeVBUMWNYREl3TWpSY1hDbGNibHh1WkdaZlFrTlRYM0J5YjNCdmNuUnBiMjV6UEMwZ1pHRjBZVjh5TURJeklDVStKVnh1SUNCdGRYUmhkR1VvWTI5dVpHbDBhVzl1WDNOamIzSmxYMjUxYldWeWFXTWdQU0JoY3k1dWRXMWxjbWxqS0dOdmJtUnBkR2x2Ymw5elkyOXlaU2tzWTI5dVpHbDBhVzl1WDNOamIzSmxYMkpwYmlBOUlHTmhjMlZmZDJobGJpaGpiMjVrYVhScGIyNWZjMk52Y21VZ1BUMGdNU0IrSUZ4Y01WeGNMQ0JqYjI1a2FYUnBiMjVmYzJOdmNtVWdQVDBnTWlCK1hGd3lYRndzWTI5dVpHbDBhVzl1WDNOamIzSmxJRDA5SURNZ2ZpQmNYRE5jWEN3Z1kyOXVaR2wwYVc5dVgzTmpiM0psSUQwOUlEUWdmbHhjTkZ4Y0xHTnZibVJwZEdsdmJsOXpZMjl5WlNBOVBTQTFJSDRnWEZ3MVhGd3NJR052Ym1ScGRHbHZibDl6WTI5eVpTQTlQU0EySUg1Y1hEWmNYQ3hqYjI1a2FYUnBiMjVmYzJOdmNtVWdQVDBnTnlCK0lGeGNOMXhjTENCamIyNWthWFJwYjI1ZmMyTnZjbVVnUFQwZ09DQitYRnc0WEZ3c0lHTnZibVJwZEdsdmJsOXpZMjl5WlNBOVBUa2dmbHhjT1Z4Y0xDQlVVbFZGSUg0Z1lYTXVZMmhoY21GamRHVnlLR052Ym1ScGRHbHZibDl6WTI5eVpTa3BLU0FsUGlWY2JpQWdaM0p2ZFhCZllua29jMmwwWlN3Z2JXOXVkR2dzSUdOdmJtUnBkR2x2Ymw5elkyOXlaVjlpYVc0cElDVStKU0JrY0d4NWNpQTZPbk4xYlcxaGNtbHpaU2hEYjNWdWREMGdiaWdwS1NBbFBpVmNiaUFnZFc1bmNtOTFjQ2dwSUNVK0pWeHVJQ0J0ZFhSaGRHVW9VSEp2Y0c5eWRHbHZiaUE5SUVOdmRXNTBMM04xYlNoRGIzVnVkQ2twWEc1Y2JtUm1YMEpEVTE5d2NtOXdiM0owYVc5dWN6d3RJRzVoTG05dGFYUW9aR1pmUWtOVFgzQnliM0J2Y25ScGIyNXpLVnh1WEc1a1psOUNRMU5mY0hKdmNHOXlkR2x2Ym5Na2JXOXVkR2dnUEMwZ1ptRmpkRzl5S0dSbVgwSkRVMTl3Y205d2IzSjBhVzl1Y3lSdGIyNTBhQ3dnYkdWMlpXeHpJRDBnWXloY1hEVmNYQ3hjWERaY1hDeGNYRGRjWEN3Z1hGdzRYRndzSUZ4Y09WeGNMQ0JjWERFd1hGd3NJRnhjTVRGY1hDa3NYRzRnSUNBZ0lDQWdJR3hoWW1Wc2N6MWpLRnhjVFdGNVhGd3NJRnhjU25WdVpWeGNMQ0JjWEVwMWJIbGNYQ3dnWEZ4QmRXZGNYQ3dnWEZ4VFpYQjBYRndzSUZ4Y1QyTjBYRndzSUZ4Y1RtOTJYRndwS1Z4dVpHWmZRa05UWDNCeWIzQnZjblJwYjI1ekpITnBkR1VnUEMwZ1ptRmpkRzl5S0dSbVgwSkRVMTl3Y205d2IzSjBhVzl1Y3lSemFYUmxMQ0JzWlhabGJITWdQU0JqS0Z4Y1FWTklRMXhjTEZ4Y1JrVk9RMXhjTEZ4Y1IwOU1SRnhjTENCY1hFeEJWVkpjWENrc1hHNGdJQ0FnSUNBZ0lHeGhZbVZzY3oxaktGeGNRWE5vSUVOeVpXVnJYRndzSUZ4Y1JtVnVZMlVnUTNKbFpXdGNYQ3dnWEZ4SGIyeGtJRk4wWVhKY1hDd2dYRnhNWVhWeVpXd2dTRzlzYkc5M1hGd3BLVnh1WEc1aVkzTXVjSEp2Y0M0eU1ESXpJRHd0SUdkbmNHeHZkQ2hrWVhSaFBXUm1YMEpEVTE5d2NtOXdiM0owYVc5dWN5d2dZV1Z6S0hnOWJXOXVkR2dzSUhrOUlGQnliM0J2Y25ScGIyNHNJR1pwYkd3OVkyOXVaR2wwYVc5dVgzTmpiM0psWDJKcGJpa3BJQ3RjYmlBZ1oyVnZiVjlpWVhJb2QybGtkR2dnUFNBdU5Td2djM1JoZEQxY1hHbGtaVzUwYVhSNVhGd3NJSEJ2YzJsMGFXOXVJRDBnWEZ4bWFXeHNYRndzSUdOdmJHOTFjaUE5SUZ4Y1lteGhZMnRjWENrcklDQmNiaUFnSUhSb1pXMWxYMkozS0NrZ0t5QWdkR2hsYldVb2NHRnVaV3d1WjNKcFpDNXRhVzV2Y2lBOUlHVnNaVzFsYm5SZllteGhibXNvS1NrcklGeHVJQ0IwYUdWdFpTaGhlR2x6TG5SbGVIUXVlQ0E5SUdWc1pXMWxiblJmZEdWNGRDaGhibWRzWlNBOUlEUTFMQ0IyYW5WemRDQTlJREVzSUdocWRYTjBQVEVzSUhOcGVtVWdQU0J5Wld3b0xqa3BLU2tyWEc0Z0lIUm9aVzFsS0hCaGJtVnNMbWR5YVdRdWJXRnFiM0l1ZUNBOVpXeGxiV1Z1ZEY5aWJHRnVheWdwS1N0Y2JpQWdiR0ZpY3loMGFYUnNaVDFjWEVKdlpIa2dZMjl1WkdsMGFXOXVJSE5qYjNKbGN5QnBiaUF5TURJelhGd3NJSGdnUFZ4Y1RXOXVkR2hjWEN3Z2VTQTlJRnhjVUhKdmNHOXlkR2x2YmlCdlppQmliMlI1SUdOdmJtUnBkR2x2YmlCelkyOXlaWE1nWEZ3cEt5QjBhR1Z0WlNoaGVHbHpMblJwZEd4bExua2dQU0JsYkdWdFpXNTBYM1JsZUhRb2MybDZaU0E5SUhKbGJDZ3hMaklwTENCaGJtZHNaU0E5T1RBcExDQmhlR2x6TG5ScGRHeGxMbmdnUFNCbGJHVnRaVzUwWDNSbGVIUW9jMmw2WlNBOUlISmxiQ2d4TGpJcExDQmhibWRzWlNBOUlEQXBLU3NnYzJOaGJHVmZabWxzYkY5aWNtVjNaWElvY0dGc1pYUjBaU0E5SUZ4Y1FteDFaWE5jWENrZ0sxeHVJQ0IwYUdWdFpTaHpkSEpwY0M1MFpYaDBJRDBnWld4bGJXVnVkRjkwWlhoMEtITnBlbVVnUFNBeE5Ta3BLMXh1SUNCMGFHVnRaU0FvYkdWblpXNWtMblJwZEd4bElEMGdaV3hsYldWdWRGOTBaWGgwS0hOcGVtVWdQU0F4T0NrcEsxeHVJQ0IwYUdWdFpTaDBhWFJzWlNBOUlHVnNaVzFsYm5SZmRHVjRkQ2h6YVhwbElEMGdNVGNwS1N0Y2JpQWdkR2hsYldVb1lYaHBjeTUwWlhoMFBXVnNaVzFsYm5SZmRHVjRkQ2h6YVhwbFBURTNLU2tyWEc0Z0lHeGhZbk1vWm1sc2JDQTlYRnhDYjJSNUlFTnZibVJwZEdsdmJpQlRZMjl5WlZ4Y0tTdGNiaUFnWm1GalpYUmZkM0poY0NoK0lITnBkR1VnS1Z4dVltTnpMbkJ5YjNBdU1qQXlNMXh1WEc0amNHUm1LSEJoYzNSbE1DaHdZWFJvSUQwZ1hGeE1ZV0pmUkdGMFlWOVVhWE56ZFdWUWNtOWpaWE56YVc1bkwyOTFkSEIxZEZ4Y0lDeGNYQzlDUTFOZmNISnZjRzl5ZEdsdmJsOHlNREl6TG5Ca1pseGNLU3dnYUdWcFoyaDBJRDBnTnl3Z2QybGtkR2dnUFNBeE15bGNiaU53Y21sdWRDaGlZM011Y0hKdmNDNHlNREl6S1Z4dUkyUmxkaTV2Wm1Zb0tWeHVYRzRqUWtOVFgzQnliM0J2Y25ScGIyNWZZV3hzSzNOallXeGxYMlpwYkd4ZmJXRnVkV0ZzS0haaGJIVmxjeUE5SUdNb1hGd2pNREF6UXpNd1hGd3NYRndqTURFMk5qVkZYRndzWEZ3ak9EQkRSRU14WEZ3c1hGd2pRemRGUVVVMVhGd3NYRndqUmpaRk9FTXpYRndzWEZ3alJFWkRNamRFWEZ3c1hGd2pRa1k0TVRKRVhGd3NJRnhjSXpoRE5URXdRVnhjTEZ4Y0l6VTBNekF3TlZ4Y0lDa3BYRzRqYzJOaGJHVmZabWxzYkY5dFlXNTFZV3dvZG1Gc2RXVnpJRDBnWXloY1hDTXdNRE5ETXpCY1hDeGNYQ013TVRZMk5VVmNYQ3hjWENNNE1FTkVRekZjWEN4Y1hDTkROMFZCUlRWY1hDeGNYR0Y2ZFhKbE1WeGNMRnhjYzJ4aGRHVm5jbUY1TVZ4Y0xGeGNjMnhoZEdWbmNtRjVNMXhjTENCY1hITnNZWFJsWjNKaGVUUmNYQ3dnWEZ4bmNtRjVNalZjWENBcEtTQmNibHh1WmlBOExTQm1kVzVqZEdsdmJpaHdZV3dwSUdKeVpYZGxjaTV3WVd3b1luSmxkMlZ5TG5CaGJDNXBibVp2VzNCaGJDd2dYRnh0WVhoamIyeHZjbk5jWEYwc0lIQmhiQ2xjYmloamIyeHpJRHd0SUdZb1hGeFpiRWR1UW5WY1hDa3BYRzVnWUdCY2JtQmdZQ0o5IC0tPlxuXG5gYGByXG5gYGByXG4jQ29kaW5nIGdyYXBocyBmb3IgTElTIENvbmZlcmVuY2UgLSBvbmx5IDIwMjMgZGF0YVxuZGF0YV8yMDIzIDwtIGRhdGFfYWxsJT4lICBmaWx0ZXIoIXllYXI9PVxcMjAyNFxcKVxuXG5kZl9CQ1NfcHJvcG9ydGlvbnM8LSBkYXRhXzIwMjMgJT4lXG4gIG11dGF0ZShjb25kaXRpb25fc2NvcmVfbnVtZXJpYyA9IGFzLm51bWVyaWMoY29uZGl0aW9uX3Njb3JlKSxjb25kaXRpb25fc2NvcmVfYmluID0gY2FzZV93aGVuKGNvbmRpdGlvbl9zY29yZSA9PSAxIH4gXFwxXFwsIGNvbmRpdGlvbl9zY29yZSA9PSAyIH5cXDJcXCxjb25kaXRpb25fc2NvcmUgPT0gMyB+IFxcM1xcLCBjb25kaXRpb25fc2NvcmUgPT0gNCB+XFw0XFwsY29uZGl0aW9uX3Njb3JlID09IDUgfiBcXDVcXCwgY29uZGl0aW9uX3Njb3JlID09IDYgflxcNlxcLGNvbmRpdGlvbl9zY29yZSA9PSA3IH4gXFw3XFwsIGNvbmRpdGlvbl9zY29yZSA9PSA4IH5cXDhcXCwgY29uZGl0aW9uX3Njb3JlID09OSB+XFw5XFwsIFRSVUUgfiBhcy5jaGFyYWN0ZXIoY29uZGl0aW9uX3Njb3JlKSkpICU+JVxuICBncm91cF9ieShzaXRlLCBtb250aCwgY29uZGl0aW9uX3Njb3JlX2JpbikgJT4lIGRwbHlyIDo6c3VtbWFyaXNlKENvdW50PSBuKCkpICU+JVxuICB1bmdyb3VwKCkgJT4lXG4gIG11dGF0ZShQcm9wb3J0aW9uID0gQ291bnQvc3VtKENvdW50KSlcblxuZGZfQkNTX3Byb3BvcnRpb25zPC0gbmEub21pdChkZl9CQ1NfcHJvcG9ydGlvbnMpXG5cbmRmX0JDU19wcm9wb3J0aW9ucyRtb250aCA8LSBmYWN0b3IoZGZfQkNTX3Byb3BvcnRpb25zJG1vbnRoLCBsZXZlbHMgPSBjKFxcNVxcLFxcNlxcLFxcN1xcLCBcXDhcXCwgXFw5XFwsIFxcMTBcXCwgXFwxMVxcKSxcbiAgICAgICAgbGFiZWxzPWMoXFxNYXlcXCwgXFxKdW5lXFwsIFxcSnVseVxcLCBcXEF1Z1xcLCBcXFNlcHRcXCwgXFxPY3RcXCwgXFxOb3ZcXCkpXG5kZl9CQ1NfcHJvcG9ydGlvbnMkc2l0ZSA8LSBmYWN0b3IoZGZfQkNTX3Byb3BvcnRpb25zJHNpdGUsIGxldmVscyA9IGMoXFxBU0hDXFwsXFxGRU5DXFwsXFxHT0xEXFwsIFxcTEFVUlxcKSxcbiAgICAgICAgbGFiZWxzPWMoXFxBc2ggQ3JlZWtcXCwgXFxGZW5jZSBDcmVla1xcLCBcXEdvbGQgU3RhclxcLCBcXExhdXJlbCBIb2xsb3dcXCkpXG5cbmJjcy5wcm9wLjIwMjMgPC0gZ2dwbG90KGRhdGE9ZGZfQkNTX3Byb3BvcnRpb25zLCBhZXMoeD1tb250aCwgeT0gUHJvcG9ydGlvbiwgZmlsbD1jb25kaXRpb25fc2NvcmVfYmluKSkgK1xuICBnZW9tX2Jhcih3aWR0aCA9IC41LCBzdGF0PVxcaWRlbnRpdHlcXCwgcG9zaXRpb24gPSBcXGZpbGxcXCwgY29sb3VyID0gXFxibGFja1xcKSsgIFxuICAgdGhlbWVfYncoKSArICB0aGVtZShwYW5lbC5ncmlkLm1pbm9yID0gZWxlbWVudF9ibGFuaygpKSsgXG4gIHRoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gNDUsIHZqdXN0ID0gMSwgaGp1c3Q9MSwgc2l6ZSA9IHJlbCguOSkpKStcbiAgdGhlbWUocGFuZWwuZ3JpZC5tYWpvci54ID1lbGVtZW50X2JsYW5rKCkpK1xuICBsYWJzKHRpdGxlPVxcQm9keSBjb25kaXRpb24gc2NvcmVzIGluIDIwMjNcXCwgeCA9XFxNb250aFxcLCB5ID0gXFxQcm9wb3J0aW9uIG9mIGJvZHkgY29uZGl0aW9uIHNjb3JlcyBcXCkrIHRoZW1lKGF4aXMudGl0bGUueSA9IGVsZW1lbnRfdGV4dChzaXplID0gcmVsKDEuMiksIGFuZ2xlID05MCksIGF4aXMudGl0bGUueCA9IGVsZW1lbnRfdGV4dChzaXplID0gcmVsKDEuMiksIGFuZ2xlID0gMCkpKyBzY2FsZV9maWxsX2JyZXdlcihwYWxldHRlID0gXFxCbHVlc1xcKSArXG4gIHRoZW1lKHN0cmlwLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSkrXG4gIHRoZW1lIChsZWdlbmQudGl0bGUgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE4KSkrXG4gIHRoZW1lKHRpdGxlID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNykpK1xuICB0aGVtZShheGlzLnRleHQ9ZWxlbWVudF90ZXh0KHNpemU9MTcpKStcbiAgbGFicyhmaWxsID1cXEJvZHkgQ29uZGl0aW9uIFNjb3JlXFwpK1xuICBmYWNldF93cmFwKH4gc2l0ZSApXG5iY3MucHJvcC4yMDIzXG5cbiNwZGYocGFzdGUwKHBhdGggPSBcXExhYl9EYXRhX1Rpc3N1ZVByb2Nlc3Npbmcvb3V0cHV0XFwgLFxcL0JDU19wcm9wb3J0aW9uXzIwMjMucGRmXFwpLCBoZWlnaHQgPSA3LCB3aWR0aCA9IDEzKVxuI3ByaW50KGJjcy5wcm9wLjIwMjMpXG4jZGV2Lm9mZigpXG5cbiNCQ1NfcHJvcG9ydGlvbl9hbGwrc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gYyhcXCMwMDNDMzBcXCxcXCMwMTY2NUVcXCxcXCM4MENEQzFcXCxcXCNDN0VBRTVcXCxcXCNGNkU4QzNcXCxcXCNERkMyN0RcXCxcXCNCRjgxMkRcXCwgXFwjOEM1MTBBXFwsXFwjNTQzMDA1XFwgKSlcbiNzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSBjKFxcIzAwM0MzMFxcLFxcIzAxNjY1RVxcLFxcIzgwQ0RDMVxcLFxcI0M3RUFFNVxcLFxcYXp1cmUxXFwsXFxzbGF0ZWdyYXkxXFwsXFxzbGF0ZWdyYXkzXFwsIFxcc2xhdGVncmF5NFxcLCBcXGdyYXkyNVxcICkpIFxuXG5mIDwtIGZ1bmN0aW9uKHBhbCkgYnJld2VyLnBhbChicmV3ZXIucGFsLmluZm9bcGFsLCBcXG1heGNvbG9yc1xcXSwgcGFsKVxuKGNvbHMgPC0gZihcXFlsR25CdVxcKSlcbmBgYFxuYGBgXG5cbjwhLS0gcm5iLXNvdXJjZS1lbmQgLS0+XG4ifQ== -->
<!-- rnb-source-begin eyJkYXRhIjoiYGBgclxuYGBgclxuI0NvZGluZyBncmFwaHMgZm9yIExJUyBDb25mZXJlbmNlIC0gb25seSAyMDIzIGRhdGFcbmRhdGFfMjAyMyA8LSBkYXRhX2FsbCU+JSAgZmlsdGVyKCF5ZWFyPT1cXDIwMjRcXClcblxuZGZfQkNTX3Byb3BvcnRpb25zPC0gZGF0YV8yMDIzICU+JVxuICBtdXRhdGUoY29uZGl0aW9uX3Njb3JlX251bWVyaWMgPSBhcy5udW1lcmljKGNvbmRpdGlvbl9zY29yZSksY29uZGl0aW9uX3Njb3JlX2JpbiA9IGNhc2Vfd2hlbihjb25kaXRpb25fc2NvcmUgPT0gMSB+IFxcMVxcLCBjb25kaXRpb25fc2NvcmUgPT0gMiB+XFwyXFwsY29uZGl0aW9uX3Njb3JlID09IDMgfiBcXDNcXCwgY29uZGl0aW9uX3Njb3JlID09IDQgflxcNFxcLGNvbmRpdGlvbl9zY29yZSA9PSA1IH4gXFw1XFwsIGNvbmRpdGlvbl9zY29yZSA9PSA2IH5cXDZcXCxjb25kaXRpb25fc2NvcmUgPT0gNyB+IFxcN1xcLCBjb25kaXRpb25fc2NvcmUgPT0gOCB+XFw4XFwsIGNvbmRpdGlvbl9zY29yZSA9PTkgflxcOVxcLCBUUlVFIH4gYXMuY2hhcmFjdGVyKGNvbmRpdGlvbl9zY29yZSkpKSAlPiVcbiAgZ3JvdXBfYnkoc2l0ZSwgbW9udGgsIGNvbmRpdGlvbl9zY29yZV9iaW4pICU+JSBkcGx5ciA6OnN1bW1hcmlzZShDb3VudD0gbigpKSAlPiVcbiAgdW5ncm91cCgpICU+JVxuICBtdXRhdGUoUHJvcG9ydGlvbiA9IENvdW50L3N1bShDb3VudCkpXG5cbmRmX0JDU19wcm9wb3J0aW9uczwtIG5hLm9taXQoZGZfQkNTX3Byb3BvcnRpb25zKVxuXG5kZl9CQ1NfcHJvcG9ydGlvbnMkbW9udGggPC0gZmFjdG9yKGRmX0JDU19wcm9wb3J0aW9ucyRtb250aCwgbGV2ZWxzID0gYyhcXDVcXCxcXDZcXCxcXDdcXCwgXFw4XFwsIFxcOVxcLCBcXDEwXFwsIFxcMTFcXCksXG4gICAgICAgIGxhYmVscz1jKFxcTWF5XFwsIFxcSnVuZVxcLCBcXEp1bHlcXCwgXFxBdWdcXCwgXFxTZXB0XFwsIFxcT2N0XFwsIFxcTm92XFwpKVxuZGZfQkNTX3Byb3BvcnRpb25zJHNpdGUgPC0gZmFjdG9yKGRmX0JDU19wcm9wb3J0aW9ucyRzaXRlLCBsZXZlbHMgPSBjKFxcQVNIQ1xcLFxcRkVOQ1xcLFxcR09MRFxcLCBcXExBVVJcXCksXG4gICAgICAgIGxhYmVscz1jKFxcQXNoIENyZWVrXFwsIFxcRmVuY2UgQ3JlZWtcXCwgXFxHb2xkIFN0YXJcXCwgXFxMYXVyZWwgSG9sbG93XFwpKVxuXG5iY3MucHJvcC4yMDIzIDwtIGdncGxvdChkYXRhPWRmX0JDU19wcm9wb3J0aW9ucywgYWVzKHg9bW9udGgsIHk9IFByb3BvcnRpb24sIGZpbGw9Y29uZGl0aW9uX3Njb3JlX2JpbikpICtcbiAgZ2VvbV9iYXIod2lkdGggPSAuNSwgc3RhdD1cXGlkZW50aXR5XFwsIHBvc2l0aW9uID0gXFxmaWxsXFwsIGNvbG91ciA9IFxcYmxhY2tcXCkrICBcbiAgIHRoZW1lX2J3KCkgKyAgdGhlbWUocGFuZWwuZ3JpZC5taW5vciA9IGVsZW1lbnRfYmxhbmsoKSkrIFxuICB0aGVtZShheGlzLnRleHQueCA9IGVsZW1lbnRfdGV4dChhbmdsZSA9IDQ1LCB2anVzdCA9IDEsIGhqdXN0PTEsIHNpemUgPSByZWwoLjkpKSkrXG4gIHRoZW1lKHBhbmVsLmdyaWQubWFqb3IueCA9ZWxlbWVudF9ibGFuaygpKStcbiAgbGFicyh0aXRsZT1cXEJvZHkgY29uZGl0aW9uIHNjb3JlcyBpbiAyMDIzXFwsIHggPVxcTW9udGhcXCwgeSA9IFxcUHJvcG9ydGlvbiBvZiBib2R5IGNvbmRpdGlvbiBzY29yZXMgXFwpKyB0aGVtZShheGlzLnRpdGxlLnkgPSBlbGVtZW50X3RleHQoc2l6ZSA9IHJlbCgxLjIpLCBhbmdsZSA9OTApLCBheGlzLnRpdGxlLnggPSBlbGVtZW50X3RleHQoc2l6ZSA9IHJlbCgxLjIpLCBhbmdsZSA9IDApKSsgc2NhbGVfZmlsbF9icmV3ZXIocGFsZXR0ZSA9IFxcQmx1ZXNcXCkgK1xuICB0aGVtZShzdHJpcC50ZXh0ID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNSkpK1xuICB0aGVtZSAobGVnZW5kLnRpdGxlID0gZWxlbWVudF90ZXh0KHNpemUgPSAxOCkpK1xuICB0aGVtZSh0aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplID0gMTcpKStcbiAgdGhlbWUoYXhpcy50ZXh0PWVsZW1lbnRfdGV4dChzaXplPTE3KSkrXG4gIGxhYnMoZmlsbCA9XFxCb2R5IENvbmRpdGlvbiBTY29yZVxcKStcbiAgZmFjZXRfd3JhcCh+IHNpdGUgKVxuYmNzLnByb3AuMjAyM1xuXG4jcGRmKHBhc3RlMChwYXRoID0gXFxMYWJfRGF0YV9UaXNzdWVQcm9jZXNzaW5nL291dHB1dFxcICxcXC9CQ1NfcHJvcG9ydGlvbl8yMDIzLnBkZlxcKSwgaGVpZ2h0ID0gNywgd2lkdGggPSAxMylcbiNwcmludChiY3MucHJvcC4yMDIzKVxuI2Rldi5vZmYoKVxuXG4jQkNTX3Byb3BvcnRpb25fYWxsK3NjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcyA9IGMoXFwjMDAzQzMwXFwsXFwjMDE2NjVFXFwsXFwjODBDREMxXFwsXFwjQzdFQUU1XFwsXFwjRjZFOEMzXFwsXFwjREZDMjdEXFwsXFwjQkY4MTJEXFwsIFxcIzhDNTEwQVxcLFxcIzU0MzAwNVxcICkpXG4jc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gYyhcXCMwMDNDMzBcXCxcXCMwMTY2NUVcXCxcXCM4MENEQzFcXCxcXCNDN0VBRTVcXCxcXGF6dXJlMVxcLFxcc2xhdGVncmF5MVxcLFxcc2xhdGVncmF5M1xcLCBcXHNsYXRlZ3JheTRcXCwgXFxncmF5MjVcXCApKSBcblxuZiA8LSBmdW5jdGlvbihwYWwpIGJyZXdlci5wYWwoYnJld2VyLnBhbC5pbmZvW3BhbCwgXFxtYXhjb2xvcnNcXF0sIHBhbClcbihjb2xzIDwtIGYoXFxZbEduQnVcXCkpXG5gYGBcbmBgYCJ9 -->
```r
```r
#Coding graphs for LIS Conference - only 2023 data
data_2023 <- data_all%>% filter(!year==\2024\)
df_BCS_proportions<- data_2023 %>%
mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ \1\, condition_score == 2 ~\2\,condition_score == 3 ~ \3\, condition_score == 4 ~\4\,condition_score == 5 ~ \5\, condition_score == 6 ~\6\,condition_score == 7 ~ \7\, condition_score == 8 ~\8\, condition_score ==9 ~\9\, TRUE ~ as.character(condition_score))) %>%
group_by(site, month, condition_score_bin) %>% dplyr ::summarise(Count= n()) %>%
ungroup() %>%
mutate(Proportion = Count/sum(Count))
df_BCS_proportions<- na.omit(df_BCS_proportions)
df_BCS_proportions$month <- factor(df_BCS_proportions$month, levels = c(\5\,\6\,\7\, \8\, \9\, \10\, \11\),
labels=c(\May\, \June\, \July\, \Aug\, \Sept\, \Oct\, \Nov\))
df_BCS_proportions$site <- factor(df_BCS_proportions$site, levels = c(\ASHC\,\FENC\,\GOLD\, \LAUR\),
labels=c(\Ash Creek\, \Fence Creek\, \Gold Star\, \Laurel Hollow\))
bcs.prop.2023 <- ggplot(data=df_BCS_proportions, aes(x=month, y= Proportion, fill=condition_score_bin)) +
geom_bar(width = .5, stat=\identity\, position = \fill\, colour = \black\)+
theme_bw() + theme(panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1, size = rel(.9)))+
theme(panel.grid.major.x =element_blank())+
labs(title=\Body condition scores in 2023\, x =\Month\, y = \Proportion of body condition scores \)+ theme(axis.title.y = element_text(size = rel(1.2), angle =90), axis.title.x = element_text(size = rel(1.2), angle = 0))+ scale_fill_brewer(palette = \Blues\) +
theme(strip.text = element_text(size = 15))+
theme (legend.title = element_text(size = 18))+
theme(title = element_text(size = 17))+
theme(axis.text=element_text(size=17))+
labs(fill =\Body Condition Score\)+
facet_wrap(~ site )
bcs.prop.2023
#pdf(paste0(path = \Lab_Data_TissueProcessing/output\ ,\/BCS_proportion_2023.pdf\), height = 7, width = 13)
#print(bcs.prop.2023)
#dev.off()
#BCS_proportion_all+scale_fill_manual(values = c(\#003C30\,\#01665E\,\#80CDC1\,\#C7EAE5\,\#F6E8C3\,\#DFC27D\,\#BF812D\, \#8C510A\,\#543005\ ))
#scale_fill_manual(values = c(\#003C30\,\#01665E\,\#80CDC1\,\#C7EAE5\,\azure1\,\slategray1\,\slategray3\, \slategray4\, \gray25\ ))
f <- function(pal) brewer.pal(brewer.pal.info[pal, \maxcolors\], pal)
(cols <- f(\YlGnBu\))
```
```
<!-- rnb-source-end -->
<!-- rnb-output-end -->
<!-- rnb-chunk-end -->
<!-- rnb-text-begin -->
<!-- rnb-text-end -->
<!-- rnb-chunk-begin -->
<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVJME52WkdsdVp5Qm5jbUZ3YUhNZ1ptOXlJRXhKVXlCRGIyNW1aWEpsYm1ObElDMGdiMjVzZVNBeU1ESXpJR1JoZEdFZ0tHTnZiblJwYm5WbFpDbGNibHh1WkdaZk0xOXdjbTl3UEMwZ2MyVnNaV04wS0dSaGRHRmZNakF5TXl3Z0xXTW9YQ0p6YVhSbFhDSXBLVnh1WEc0alEyOTFiblFnWVc1a0lIQnliM0J2Y25ScGIyNGdkR0ZpYkdVZ1ptOXlJREl3TWpNZ1ltTnpJSE5qYjNKbGN5QmllU0J0YjI1MGFGeHVZbU56WDNSaFlteGxJRHd0SUhSaFlteGxLR1JoZEdGZk1qQXlNeVJqYjI1a2FYUnBiMjVmYzJOdmNtVXNJR1JoZEdGZk1qQXlNeVJ0YjI1MGFDd2daR0YwWVY4eU1ESXpKSE5wZEdVcFhHNXVZVzFsY3loa2FXMXVZVzFsY3loaVkzTmZkR0ZpYkdVcEtTQThMU0JqS0Z3aVFtOWtlU0JEYjI1a2FYUnBiMjRnVTJOdmNtVmNJaXdnWENKTmIyNTBhRndpTENCY0lsTnBkR1ZjSWlsY2JtSmpjMTkwWVdKc1pTQThMU0JoWkdSdFlYSm5hVzV6S0dKamMxOTBZV0pzWlNsY2JtSmpjMTkwWVdKc1pWeHVYRzVnWUdBaWZRPT0gLS0+XG5cbmBgYHJcbiNDb2RpbmcgZ3JhcGhzIGZvciBMSVMgQ29uZmVyZW5jZSAtIG9ubHkgMjAyMyBkYXRhIChjb250aW51ZWQpXG5cbmRmXzNfcHJvcDwtIHNlbGVjdChkYXRhXzIwMjMsIC1jKFxcc2l0ZVxcKSlcblxuI0NvdW50IGFuZCBwcm9wb3J0aW9uIHRhYmxlIGZvciAyMDIzIGJjcyBzY29yZXMgYnkgbW9udGhcbmJjc190YWJsZSA8LSB0YWJsZShkYXRhXzIwMjMkY29uZGl0aW9uX3Njb3JlLCBkYXRhXzIwMjMkbW9udGgsIGRhdGFfMjAyMyRzaXRlKVxubmFtZXMoZGltbmFtZXMoYmNzX3RhYmxlKSkgPC0gYyhcXEJvZHkgQ29uZGl0aW9uIFNjb3JlXFwsIFxcTW9udGhcXCwgXFxTaXRlXFwpXG5iY3NfdGFibGUgPC0gYWRkbWFyZ2lucyhiY3NfdGFibGUpXG5iY3NfdGFibGVcblxuYGBgXG5cbjwhLS0gcm5iLXNvdXJjZS1lbmQgLS0+XG4ifQ== -->
#Coding graphs for LIS Conference - only 2023 data (continued)
df_3_prop<- select(data_2023, -c(\site\))
#Count and proportion table for 2023 bcs scores by month
bcs_table <- table(data_2023$condition_score, data_2023$month, data_2023$site)
names(dimnames(bcs_table)) <- c(\Body Condition Score\, \Month\, \Site\)
bcs_table <- addmargins(bcs_table)
bcs_table
<!-- rnb-output-end -->
<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVZR0JnY2x4dUkwTnZaR2x1WnlCbmNtRndhSE1nWm05eUlFeEpVeUJEYjI1bVpYSmxibU5sSUMwZ2IyNXNlU0F5TURJeklHUmhkR0VnS0dOdmJuUnBiblZsWkNsY2JseHVaR1pmTTE5d2NtOXdQQzBnYzJWc1pXTjBLR1JoZEdGZk1qQXlNeXdnTFdNb1hGeHphWFJsWEZ3cEtWeHVYRzRqUTI5MWJuUWdZVzVrSUhCeWIzQnZjblJwYjI0Z2RHRmliR1VnWm05eUlESXdNak1nWW1OeklITmpiM0psY3lCaWVTQnRiMjUwYUZ4dVltTnpYM1JoWW14bElEd3RJSFJoWW14bEtHUmhkR0ZmTWpBeU15UmpiMjVrYVhScGIyNWZjMk52Y21Vc0lHUmhkR0ZmTWpBeU15UnRiMjUwYUN3Z1pHRjBZVjh5TURJekpITnBkR1VwWEc1dVlXMWxjeWhrYVcxdVlXMWxjeWhpWTNOZmRHRmliR1VwS1NBOExTQmpLRnhjUW05a2VTQkRiMjVrYVhScGIyNGdVMk52Y21WY1hDd2dYRnhOYjI1MGFGeGNMQ0JjWEZOcGRHVmNYQ2xjYm1KamMxOTBZV0pzWlNBOExTQmhaR1J0WVhKbmFXNXpLR0pqYzE5MFlXSnNaU2xjYm1KamMxOTBZV0pzWlZ4dVhHNWdZR0JjYm1CZ1lDSjkgLS0+XG5cbmBgYHJcbmBgYHJcbiNDb2RpbmcgZ3JhcGhzIGZvciBMSVMgQ29uZmVyZW5jZSAtIG9ubHkgMjAyMyBkYXRhIChjb250aW51ZWQpXG5cbmRmXzNfcHJvcDwtIHNlbGVjdChkYXRhXzIwMjMsIC1jKFxcc2l0ZVxcKSlcblxuI0NvdW50IGFuZCBwcm9wb3J0aW9uIHRhYmxlIGZvciAyMDIzIGJjcyBzY29yZXMgYnkgbW9udGhcbmJjc190YWJsZSA8LSB0YWJsZShkYXRhXzIwMjMkY29uZGl0aW9uX3Njb3JlLCBkYXRhXzIwMjMkbW9udGgsIGRhdGFfMjAyMyRzaXRlKVxubmFtZXMoZGltbmFtZXMoYmNzX3RhYmxlKSkgPC0gYyhcXEJvZHkgQ29uZGl0aW9uIFNjb3JlXFwsIFxcTW9udGhcXCwgXFxTaXRlXFwpXG5iY3NfdGFibGUgPC0gYWRkbWFyZ2lucyhiY3NfdGFibGUpXG5iY3NfdGFibGVcblxuYGBgXG5gYGBcblxuPCEtLSBybmItc291cmNlLWVuZCAtLT5cbiJ9 -->
<!-- rnb-source-begin eyJkYXRhIjoiYGBgclxuYGBgclxuI0NvZGluZyBncmFwaHMgZm9yIExJUyBDb25mZXJlbmNlIC0gb25seSAyMDIzIGRhdGEgKGNvbnRpbnVlZClcblxuZGZfM19wcm9wPC0gc2VsZWN0KGRhdGFfMjAyMywgLWMoXFxzaXRlXFwpKVxuXG4jQ291bnQgYW5kIHByb3BvcnRpb24gdGFibGUgZm9yIDIwMjMgYmNzIHNjb3JlcyBieSBtb250aFxuYmNzX3RhYmxlIDwtIHRhYmxlKGRhdGFfMjAyMyRjb25kaXRpb25fc2NvcmUsIGRhdGFfMjAyMyRtb250aCwgZGF0YV8yMDIzJHNpdGUpXG5uYW1lcyhkaW1uYW1lcyhiY3NfdGFibGUpKSA8LSBjKFxcQm9keSBDb25kaXRpb24gU2NvcmVcXCwgXFxNb250aFxcLCBcXFNpdGVcXClcbmJjc190YWJsZSA8LSBhZGRtYXJnaW5zKGJjc190YWJsZSlcbmJjc190YWJsZVxuXG5gYGBcbmBgYCJ9 -->
```r
```r
#Coding graphs for LIS Conference - only 2023 data (continued)
df_3_prop<- select(data_2023, -c(\site\))
#Count and proportion table for 2023 bcs scores by month
bcs_table <- table(data_2023$condition_score, data_2023$month, data_2023$site)
names(dimnames(bcs_table)) <- c(\Body Condition Score\, \Month\, \Site\)
bcs_table <- addmargins(bcs_table)
bcs_table
```
```
<!-- rnb-source-end -->
<!-- rnb-output-end -->
<!-- rnb-chunk-end -->
<!-- rnb-text-begin -->
<!-- rnb-text-end -->
<!-- rnb-chunk-begin -->
<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVjMlZ6YzJsdmJrbHVabThvS1Z4dVlHQmdJbjA9IC0tPlxuXG5gYGByXG5zZXNzaW9uSW5mbygpXG5gYGBcblxuPCEtLSBybmItc291cmNlLWVuZCAtLT5cbiJ9 -->
sessionInfo()
````
```r
sessionInfo()
```
Add a new chunk by clicking the Insert Chunk button on the toolbar or by pressing Ctrl+Alt+I.
When you save the notebook, an HTML file containing the code and output will be saved alongside it (click the Preview button or press Ctrl+Shift+K to preview the HTML file).
The preview shows you a rendered HTML copy of the contents of the editor. Consequently, unlike Knit, Preview does not run any R code chunks. Instead, the output of the chunk when it was last run in the editor is displayed.